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Detailed  instructions  are  presented  for  using  a  computer  program  which 
calculates  the  six-degree-of-freedom  trajectories  of  external  stores  which  are 
separated  from  fighter-bomber  type  aircraft  flying  at  supersonic  speeds. 

Multiple  circular  or  elliptic  store  configurations  may  be  handled.  Parent 
aircraft  configurations  may  consist  of  a  circular  or  arbitrary  cross  suction 
fuselage  with  ramp  external  compression  inlets,  and  a  wing,  pylon,  and  rack. 

The  program  uses  linear  potential-flow  theory  to  model  the  wing  and  pylon 
loading  and  thickness.  Three-dimensional  line  sources  and  doublets  are  used  5 
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to  model  circular  fuselages  and  stores.  The  poncircular  fuselage  and 
elliptic  store  surfaces  are  modeled  with  constant  source  panels.  Nonlinear 
corrections  are  made  to  the  wing,  fuselage,  rack,  store  and  fuselage  inlet 
models  to  simulate  shocks.  The  program  alsp  calculates  the  trajectory  of 
the  store  as  it  separates  from  the  aircraft.  This  report  describes  the 
program,  presents  instructions  for  preparing  input  for  the  program,  describes 
the  output  from  the  program,  and  presents  a  sample  case.  The  program  repre¬ 
sents  an  extension  of  an  earlier  program  restricted  to  circular  bodies  at 
supersonic  speeds  written  by  the  present  authors  and  described  in  AFFDL-TR- 
76-41. 

This  volume  presents  the  detailed  descriptions  of  the  calculations 
performed  in  each  of  the  subroutines  in  Program  I.  Alsp  included  are  the 
descriptions  of  each  of  the  variables  passed  between  routines. 
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FOREWORD 


This  report  entitled  "Prediction  of  Supersonic  Store  Separation 
Characteristics  Including  Fuselage  and  Stores  of  Noncircular  Cross 
Section,"  describes  a  combined  theoretical-experimental  program 
directed  toward  developing  a  computer  program  for  predicting  the 
trajectory  of  an  external  store  separated  from  an  aircraft  flying 
at  supersonic  speed.  It  represents  an  extension  of  previous  work 
covered  in  AFFDL-TR-76-41  to  include  more  realistic  modeling  of 
fuselage  shapes  including  noncircular  cross  sections  and  ramp  type 
engine  air  inlets,  and  to  include  modeling  store  shapes  with 
elliptic  cross  section  with  multiple  sets  of  arbitrarily  oriented 
fins.  Volume  I,  "Theoretical  Methods  and  Comparisons  with  Experi¬ 
ment",  describes  the  theoretical  approach  and  presents  extensive 
comparisons  with  experimental  data.  Volume  II,  "Users  Manual  for 
the  Computer  Program" ,  presents  detailed  instructions  on  the  use 
of  the  computer  program  with  emphasis  on  preparation  of  input  data 
and  interpretation  of  output.  This  volume,  Volume  III,  "Appendices 
A  and  B,  Details  of  Program  I",  provides  additional  descriptions 
of  the  individual  subroutines  and  program  variables  passed  between 
modules  in  the  first  of  two  programs.  Volume  IV,  "Appendices  C 
and  D,  Details  of  Program  II",  provides  additional  descriptions  of 
the  individual  subroutines  and  program  variables  passed  between 
modules  in  the  second  program. 

This  work  was  carried  out  by  Nielsen  Engineering  &  Research, 
Inc.  (NEAR),  510  Clyde  Avenue,  Mountain  View,  California  94043, 
under  Contract  No.  F33615-76-C-3077.  The  contract  was  initiated 
under  Project  2403,  Task  240305,  of  the  Air  Force  Flight  Dynamics 
Laboratory.  The  Air  Force  Project  Engineer  on  the  contract  was 
Calvin  L.  Dyer,  AFWAL/FIGC.  The  report  number  assigned  by  Nielsen 
Engineering  &  Research,  Inc.  is  NEAR  TR  210. 
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PREDICTION  OF  SUPERSONIC  STORE  SEPARATION 
CHARACTERISTICS  INCLUDING  FUSELAGE  AND 
STORES  OF  NONCIRCULAR  CROSS  SECTION 

Volume  III  -  Appendices  A  and  B,  Details  of  Program  I 

SUMMARY 

The  purpose  of  this  volume  is  to  provide  the  additional 
details  of  the  parts  of  the  first  program  that  would  be  useful 
to  the  programmer  or  engineer  interested  in  understanding  the 
calculations  and  computer  code  herein.  The  information  included 
here  as  Appendices  A  and  B  describes  the  function  and  operations 
performed  by  each  routine,  a  description  of  data  transferred 
between  subroutines  and  a  listing  of  Program  I  itself. 

Appendix  A  provides  a  detailed  description  of  the  operations 
and  flow  of  calculations  of  each  of  the  individual  routines  in 
Program  I.  Included  are  a  description  of  the  flow  of  the  calcu¬ 
lations,  including  flow  charts  of  some  routines,  a  description 
of  any  program  arguments,  and  a  program  listing. 

Appendix  B  provides  a  description  of  all  variables  passed 
between  routines  in  common  blocks.  A  listing  of  each  common 
block  with  a  description  of  each  variable,  array,  or  index  in 
the  common  is  provided.  A  special  section  is  provided  for 
the  multiple  uses  of  blank  common  as  well  as  a  cross  reference 
chart  of  routine  versus  common  block  usage. 


APPENDIX  A 

DETAILS  OF  PROGRAM  I  SUBROUTINES 

A-I  Introduction 


The  purpose  of  this  appendix  is  to  provide  more  detailed 
information  on  the  calculations  performed  in  the  first  program 
which  was  described  in  Section  3  of  Volume  II.  A  listing  of 
Program  I  is  presented  in  Figure  A-l  and  a  general  flow  chart  in 
Figure  A-2.  This  appendix  will  present  the  details  of  the  flow  of 
the  calculations  for  each  routine,  a  description  of  the 
variables  in  the  argument  list,  individual  detailed  flow  charts, 
and  a  summary  of  the  functions  of  each  routine.  The  program 
consists  of  a  main  program  and  59  subroutines.  The  main  program 
will  be  described  and  then  the  subroutines  will  be  described  in 
alphabetical  order.  The  subroutines  and  their  functions  are 
listed  in  Table  A-l.  Flow  charts  of  some  of  the  individual 
routines  are  shown  in  Figures  A-3  through  A-18.  Refer  to  Volume 
II  of  this  report  for  the  list  of  symbols  used  in  this  text. 


TABLE  A-l.  SUBROUTINES  USED  IN  PROGRAM  I 


Subroutine 

Name 

Function 

LDCALC 

main  program  to  organize  calculation  flow 

ASECTN 

calculates  cross  sectional  area  of  an  arbitrarily 

paneled  body 

BDYGEN 

calculates  line  sources  and  doublets  to  give  a 
required  body  shape  and  angle  of  attack 

BLYOUT 

lays  out  fuselage  constant  u-velocity  panels  for 

circular  body 

BLY0T2 

lays  out  noncircular  fuselage  constant  u-velocity 

panels 

2 


BOD PAN 


BODPN2 

BODVEL 

BSHOCK 

CONFIG 

DOUBLT 

DPCOEF 

DPRHS 

FLDVEL 

FLDVL2 

FRSTRT 

FUSEIO 

GEOM 


organizes  the  calculation  ot  the  revised  axial 
spacing  foi  the  noncircular  fuselage  source  panels 
performs  the  revision  of  the  axial  spacing  on  the 
noncircular  bodies  and  computes  the  body  panel 
geometry 

organizes  the  looping  through  source  panels  to 

compute  panel  influence  coefficients 

computes  the  nonlinear  shock  wave  shape  emanating 

from  the  nose  of  a  noncircular  body 

reads  and  prints  the  input  geometrical  description 

of  the  external  shape  of  noncircular  fuselage  or 

elliptic  stores 

calculates  the  strength  of  a  linear  line  doublet 
calculates  coefficient  matrix  of  the  set  of  equations 
to  be  solved  for  the  wing-pylon-fuselage  constant 
u-velocity  panel  strengths 

calculates  right-hand  side  of  the  set  of  equations 
to  be  solved  for  the  wing-pylon-fuselage  constant 
u-velocity  panel  strengths 

organizes  the  calculation  of  the  u,v,w  velocity 
components  at  the  field  points  XFP,YFP,ZFP 
computes  the  three  components  of  velocity  induced 
at  field  points  by  the  source  panels  of  a  given 
body  ring 

saves  or  restores  required  program  information 
for  source  panel  model  to  restart  configuration 
analysis 

reads  and  prints  fuselage  data  and  calls  BDYGEN 
to  calculate  line  source  and  doublet  distributions 
and  SHKSHP  to  compute  nonlinear  shock  for  circular 
bod i es 

organizes  the  reading  and  printing  of  the  noncircular 
body  input  and  the  calculation  of  the  source  panel 
geometry  arrays 
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INLSHK 


INLTST 

INLXYZ 

INVER1 

I OREAD 

IOWRIT 

NETCLC 

NEW RAD 

NEWRD2 

NULYT 

PANEL 

PANVEL 

PAS001 

PAS002 

PLYOUT 


computes  the  nonlinear  shock  shape  emanating 
from  a  ramp  type  inlet 

determines  whether  a  source  panel  is  an  inlet 
panel 

scans  the  x,y,z  coordinates  of  the  inlet  panel 
corners  to  determine  inlet  geometric  parameters 
solves  a  system  of  simultaneous  linear  algebraic 
equations 

performs  an  unformatted  read  from  external 
file,  10 

performs  an  unformatted  write  onto  the  external 
file,  10 

calculates  the  net  corner  strengths  for  a  given 
set  of  wing,  pylon,  or  fuselage  constant  u-velocity 
panels 

organizes  the  revision  of  the  noncircular  body 
meridian  line  spacing  for  the  source  panels 
computes  the  revised  noncircular  body  meridional 
line  spacing  for  a  given  body  segment 
selects  a  subset  of  constant  u-velocity  and  wing 
and  pylon  thickness  panel  layout  data  and  calculates 
and  saves  the  net  strengths  associated  with  each 
panel  corner 

calculates  the  direction  cosines  of  the  normal 
vector,  the  centroid,  area  and  inclination  angles 
of  an  arbitrary  quadrilateral  panel 
transforms  field  points  into  local  panel 
coordinates  and  calls  SORPAN  for  the  calculation 
of  the  panel  influence  coefficient 
performs  the  L*U  decomposition  of  a  positive 
definite  matrix 

solves  the  system  of  equations  [L*U]*X  =  B 
reads  and  prints  pylon  data  and  lays  out 
constant  u-velocity  panels  on  the  pylon 
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PYTHIN 
RAC  K 1 0 

SHAPE 

SHKSHP 

SMARC  H 

SOLVE 

SORPAN 

SOURC E 
STORIO 

SWNGIN 

THKLYT 

THKOUT 

VELBD1 

VELCAL 


reads  pylon  thickness  data 

reads  and  prints  the  rack  data  and  calculates 
the  rack  source  and  doublet  distributions 
calculates  the  radius  and  surface  slope  at  a 
point  on  a  body  from  input  polynomials 
calculates  the  nonlinear  shape  of  the  shock  wave 
shape  produced  by  an  axisymmetric  body  at  zero 
angle  of  attack 

solves  for  the  source  panel  strengths  in 
supersonic  flow  using  a  ring  by  ring  marching 
technique 

computes  the  source  panel  boundary  condition  for 
a  noncircular  body  in  the  presence  of  a  free  stream 
angle  of  attack  and  calls  SMARCH  to  compute  the 
panel  strengths 

computes  the  three  velocity  components  induced  at 
a  specified  control  point  by  a  constant  source  panel 
inclined  to  the  free  stream 

calculates  the  strength  of  a  linear  line  source 
reads  and  prints  store  data  and  computes  circular 
body  line  source  and  doublet  distributions  and 
shock  shapes  or  elliptic  store  source  panel 
strengths  and  shock  shape 

reads  and  prints  wing  constant  u-velocity  panel 
data  and  twist  and  camber  data,  if  any 
lays  out  wing  and  pylon  thickness  panels 
prints  wing  and  pylon  thickness  panel  input  data 
calculates  influence  functions  due  to  a  fuselage 
constant  u-velocity  panel 

calculates  velocities  at  a  field  point  due  to 
fuselage  or  store  line  sources  and  doublets 
organizes  and  calls  for  the  calculation  of  the 
aerodynamic  influence  coefficient  matrices  for  the 
body  source  panels 


VELCMP 


VEL01 


VELOT1 

VELPP1 

VELPT1 

VELWP1 

VELWT1 

WDYBDY 

WITHIN 

WLYOUT 

WRFILE 

YZBIP 

YZBIP2 


calculates  the  influence  of  a  semi-infinite 
triangle  associated  with  a  constant  u-velocity 
panel 

calculates  the  influence  of  a  semi-infinite 
triangle  associated  with  a  thickness  source  panel 
calculates  influence  functions  due  to  a  pylon 
constant  u-velocity  panel 

calculates  velocities  at  a  field  point  due  to 
pylon  thickness  source  panels 
calculates  influence  functions  due  to  a  wing 
constant  u-velocity  panel 

calculates  velocities  at  a  field  point  due  to 
wing  thickness  source  panels 

reads  and  prints  noncircular  fuselage  data  and 
organizes  the  calls  for  the  body  source  panel 
layout,  source  strengths  and  nonlinear  shock 
shape 

reads  wing  thickness  data 

lays  out  constant  u-velocity  panels  on  wing 
writes  the  dataset  on  TAPE12  which  contains  the 
information  required  to  continue  the  supersonic 
store  separation  computations  in  Program  II 
organizes  the  scan  of  the  source  panel  geometry  to 
define  the  cross  section  of  the  noncircular 
fuselage  interference  shell 

interpolates  in  the  panel  geometry  of  a  segment  to 
define  the  y-z  values  of  the  noncircular  interference 
shell 
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A-2  Program  LDCALC 


1 


The  flow  chart  which  was  presented  in  Figure  2  of  Volume  II 
details  the  basic  flow  of  the  calculations  in  program  LDCALC 
and  will  not  be  expanded  in  this  appendix.  A  written  summary  of 
the  calculations  to  be  made  in  the  first  program  is  presented. 

A  listing  of  the  program  is  presented  in  Figures  A-l (a)  and  A-l (b) 
of  this  report. 

LDCALC  organizes  the  reading  and  printing  of  input  and  the 
calculation  of  the  strengths  of  the  singularities  modeling  the 
individual  parent  aircraft  and  store  components.  The  constants, 
flow  conditions,  and  component  options  are  first  defined.  The 
source  strength  solutions  are  then  defined  for  the  circular 
fuselage  in  FUSEIO  or  for  the  noncircular  fuselage  in  WDYBDY.  The 
wing  u-velocity  panel  data  are  read  in  SWNGIN  and  the  panel  layout 
is  defined  in  WLYOUT.  The  wing  thickness  data  are  read  in  WITHIN. 
PLYOUT  and  PYTHIN  perform  similar  functions  for  the  pylon.  The 
thickness  distribution  is  printed  from  THKOUT  and  the  thickness 
panel  layout  for  wing  and  pylon  occurs  in  THKLYT .  The  u-velocity 
panel  layout  for  the  circular  fuselage  is  generated  in  BLYOUT  and 
in  BLY0T2  for  the  noncircular  fuselage.  A  source  strength  model 
may  also  be  generated  for  the  rack  in  RACKIO.  The  source  strength 
determination  for  either  circular  or  elliptic  store  bodies  are 
then  computed  in  STORIO.  The  right  hand  side  and  influence 
coefficients  for  the  u-velocity  panels  are  computed  in  DPRHS  and 
DPCOEF,  respectively.  The  panel  solutions  are  solved  for  from  the 
set  of  simultaneous  equations  in  INVER1.  A  summary  of  the  panel 
control  points  and  solutions  is  printed.  The  modified  layout 
based  on  strengths  summed  at  panel  corners  is  then  generated  in 
NULYT.  The  necessary  arrays  and  variables  are  lastly  saved  on 
TAPE12  in  WRFILE  for  restarting  the  problem  in  the  second  program. 


Subroutine  references 

BLYOUT,  BLYOT2,  DPCOEF,  DPRHS ,  FUSEIO,  NULYT,  PLYOUT,  PYTHIN, 
RAC  K 1 0 ,  STORIO,  SWNGIN,  THKLYT ,  THKOUT,  WDYBDY,  WITHIN,  WLYOUT, 
WRFILE. 


A-3  Subroutine  ASECTN 


Subroutine  ASECTN  computes  the  cross  sectional  area  of  an 
arbitrarily  paneled  body.  It  is  used  to  compute  the  area  enclosed 
by  arbitrary  Y  and  Z  values  at  each  station  in  the  external 
specification  of  the  body  shape.  It  is  called  from  CONFIG. 

The  equation  used  to  define  the  area  is 


NRAD 


AREA  = 


l  4  (Y .  •  Z  .  .  -  Y.  .  •  Z 

L  _  2  1  l-l  l-l  l 


( A-l ) 


i  =  2 


If  the  body  is  symmetric  the  sum  is  computed  for  the  half  body 
and  doubled.  A  listing  of  the  routine  is  presented  in  Figure  A-l(b) 
of  this  report.  A  description  of  the  parameters  in  the  argument 
list  follows: 


arrays  of  section  corner  points 
number  of  points  around  section 
area  of  cross  section 

Called  by 
CONFIG 


Y,Z 

NRAD 

AREA 


A-4  Subroutine  BDYGEN 


Subroutine  BDYGEN  calculates  the  line  source  and  doublet 
strengths  using  control  points  on  the  surface  of  the  fuselage 
and  store.  The  method  used  is  described  in  Appendix  I  of 
Reference  1.  A  listing  of  the  subroutine  is  presented  in 
Figure  A-l (c)  and  a  flow  chart  in  Figure  A-3  of  this  report.  The 
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coordinate  system  associated  with  the  subroutine  is  shown  in  the 
sketches  of  Appendix  I  of  Reference  1. 


At  the  beginning  of  the  subroutine  a  test  is  performed  to 
determine  if,  at  the  base  of  the  body,  the  radial  distance  ti¬ 
the  Mach  cone  emanating  from  the  body  nose  is  less  than  the 
maximum  radius  of  the  body.  If  so,  an  error  message  is  printed 
out  (see  Section  3.5  of  Volume  II)  and  the  program  stops. 

Next,  as  the  flow  chart  indicates,  N  is  set  equal  to  NXBODY  -  1 
and  the  body  axis  is  divided  into  N  segments  of  equal  length. 

The  x  locations  of  the  body  definition  points,  XBODY(J),  are 
determined  at  these  equally  spaced  axis  points  and,  inside  a 
short  DO  loop,  subroutine  SHAPE  is  called  to  calculate  the 
radius  of  the  body,  RBODY,  and  the  surface  slope,  RPBODY,  at 
each  body  definition  point.  Subroutine  SHAPE  requires  that 
the  shape  definition  quantities  be  made  dimensionless  by  body 
length.  Subroutine  BDYGEN  accounts  for  this  before  and  after 
the  calls  to  subroutine  SHAPE.  Next,  the  control  points  are 
located  effectively  midway  between  the  body  definition  points. 
Subroutine  SHAPE  calculates  the  body  radius,  RF,  and  the 
surface  slope,  DRDX,  at  each  control  point.  Finally,  the  axis 
points,  TX,  which  are  the  origins  of  the  conical  line  sources 
and  doublets,  are  determined. 

The  next  sections  of  the  subroutine  are  devoted  to  revising 
the  layout  of  the  body  definition  points  and  control  points  and 
the  origins  of  the  line  singularities  if  the  first  control  point 
lies  outside  the  Mach  cone  originating  at  the  body  nose.  An 
iteration  is  performed  to  determine  the  intersection  of  the  cone 
with  the  body  surface.  Once  this  point  is  found  the  body 
definition  points  are  redistributed  over  the  remainder  of  the 
body.  The  procedure  described  in  the  previous  paragraph  is  used 
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to  redefine  the  control  points  and  the  origins  of  the  line 
singularities. 

The  remainder  of  the  subroutine  calculates  the  source  and 

doublet  strengths  at  the  control  points.  Subroutine  SOURCE  is 

called  and  the  source  strength  at  the  first  control  point  is 

calculated  using  Equation  (1-14)  and  at  the  remaining  control 

points  using  Equation  (1-17)  of  Reference  1.  One  should  note 

that  the  Ith  source  strength,  T(I),  is  the  constant  in  the 

algebraic  notation  of  the  equations.  Similarly,  using  subroutine 

DOUBLT,  the  doublet  strengths  are  calculated  using  Equation  (1-27) 

for  the  first  control  point  and  Equation  (1-28)  of  Reference  1  for 

the  remaining  control  points.  Here,  the  Ith  doublet  strength, 

TC  ( I )  ,  equals  K,,  „  .  . 

d,  1-1 

After  calculating  the  source  and  doublet  strengths,  the 
subroutine  prints  the  body  definition  point  data,  the  singularity 
origins,  and  the  singularity  strengths  and  sums  the  source  and 
doublet  strengths.  A  description  of  the  parameters  in  the 
argument  list  follow: 


NXBODY  number  of  body  definition  points 

RADIUS  maximum  radius  of  body 


BODYL 

NSEG 

XEND 


SCOEF 


length  of  actual  body 

number  of  polynomials  used  to  specify  body  shape 

array  containing  x  locations  of  end  points  of 
polynomial  sections  specifying  body  shape,  made 
dimensionless  by  body  length 

array  containing  coefficients  of  polynomials 
used  to  specify  body  shape 
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T 


TC 

TX 

ALPHA 

RPBODY(I) 

LBODY 

SUMK 

SUMKD 

Subroutine 

DOUBLT 

Called  by 
FUSE 10 


array  containing  source  strengths;  T(J)  is  the 
quantity  Kj^  calculated  by  Equations  (1-14)  and 
(1-17)  of  Reference  1 

array  containing  doublet  strengths;  TC(J)  is  the 
quantity  calculated  by  Equations  (1-27)  and 

(1-28)  of  Reference  1 

array  containing  x  locations  of  origins  of  conical 
line  sources  and  doublets;  positive,  measured  from  tip 
of  nose 

angle  measured  between  the  body  centerline  and  the 
free-stream  direction 

slope  of  body  surface  at  Ith  body  definition  point 

BODYL*XEND (NSEG) 

sum  of  the  source  strengths 

sum  of  the  doublet  strengths 

references 
SHAPE,  SOURCE 

RAC  K 1 0 ,  STORIO 


A-5  Subroutine  BLYOUT 


Subroutine  BLYOUT  calculates  quantities  needed  to  define  the 
constant  u-velocit"  panels  on  the  circular  fuselage  and  to 
locate  the  panel  control  points.  The  corner  and  control  point 
coordinates  calculated  for  the  body  panels  are  stored  in  arrays 
in  locations  which  follow  the  same  quantities  calculated  for 
the  wing  and  pylon  constant  u-velocity  panels.  All  coordinates 
are  in  the  wing  coordinate  system  which  is  shown  in  Figure  6 
of  Volume  II.  A  listing  of  the  subroutine  is  presented  in 
Figure  A-l (d)  and  a  flow  chart  in  Figure  A-4 . 

The  fuselage  interference  panels  are  laid  out  on  the 
left-half  of  the  fuselage  surface,  as  shown  in  Figure  3  of  Volume 
II.  The  panel  numbering  convention  employed  in  a  typical  ring 
containing  five  panels,  three  above  the  wing  and  two  below,  is 
shown  in  the  following  sketch. 


At  the  beginning  of  the  subroutine,  quantities  DX,  NBD,  NBIP, 
and  ANGW,  associated  with  the  geometry  of  the  entire  panel 
layout,  are  calculated.  Next,  the  variable  ANGW  is  tested  to 
determine  if  th  ving  is  tangent  to  the  fuselage,  either  at 
the  top  or  at  the  bottom.  As  the  flow  chart  shows,  this  tangency 
condition  determines  the  initialization  procedure  performed  before 
the  main  loop  calculations. 

The  main  part  of  the  subroutine  is  a  double  DO  loop  with  the 
outer  loop  index  running  over  the  rows  above  or  below  the  wing 
and  the  inner  loop  index  over  the  NCWB  panels  in  a  lengthwise 
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row.  If  the  wing  is  not  tangent,  this  double  loop  is  performed 
first  for  the  panels  below  the  wing,  then  repeated  for  the 
panels  above  the  wing. 

For  a  given  lengthwise  row  the  y,z  coordinates  of  the  panel 
corners  and  the  control  points  are  first  obtained  in  the  panel 
coordinate  system.  The  following  sketch  shows  the  coordinate 
system  and  certain  quantities  associated  with  a  panel  located 
above  the  wing. 


Within  the  inner  DO  loop  the  coordinates  are  transformed  into  the 
wing  coordinate  system  and  stored  in  appropriate  array  locations. 
Polar  angle  THTI  and  quantities  SNT2  and  CST2,  necessary  for 
panel-wing  transformations,  are  saved  in  arrays  for  subsequent 
use  by  other  routines.  The  corner  and  control  point  xw  coordinates 
are  calculated  and  stored;  no  transformations  are  necessary  for 
these  quantities. 

The  convention  employed  in  labeling  the  corner  coordinates  for 
the  fuselage  interference  panels  is  that  the  right  side  of  the 
panel  is  located  in  clockwise  rotation  from  the  left  side,  when 
looking  forward. 

Called  by 

LDCALC 


A-6  Subroutine  BLY0T2 


Subroutine  BLY0T2  computes  the  u-velocity  panel  corner  points, 
control  points,  and  transformation  angles  for  the  noncircular 
fuselage  body  shape.  The  section  geometry  is  obtained  from  YZBIP 
at  the  user  specified  body  x-station  in  the  arbitrary  fuselage 
geometry  description.  The  section  generated  is  held  constant  in 
cross  section  in  the  axial  direction.  No  logic  has  been  provided 
to  force  panel  edges  at  the  wing-body  interface  to  match.  The 
initial  paneling  of  the  body  by  the  user  must  be  performed  so  as 
to  guarantee  the  matching  of  the  panel  edges  with  the  wing  root 
chord.  A  listing  of  the  routine  is  presented  in  Figure  A-l(e)  of 
this  report. 

BLY0T2  performs  two  loops  when  generating  the  tables  of  panel 
properties.  The  outer  DO  loop  steps  through  the  number  of  panels 
around  the  circumference  to  generate  the  local  section  properties 
as  shown  in  Figure  18  of  Volume  II.  The  inner  DO  loop  increments 
the  indices  axially  to  repeat  the  identical  section  geometry  for 
the  remaining  downstream  panels  on  the  body  interference  shell. 

The  panel  properties  are  laid  out  starting  at  the  bottom  for  the 
left  hand  side.  All  quantities  defined  maintain  the  same 
definitions  and  conventions  used  in  BLYOUT  for  the  circular 
fuselage. 

Called  by 
LDCALC 

A-7  Subroutine  BODPAN 

Subroutine  BODPAN  organizes  the  calculation  of  the  revised 
axial  spacing  for  noncircular  fuselage  source  panels.  The 
starting  locations  in  blank  common  of  all  arrays  used  to  define 
the  geometric  properties  of  the  panels  are  computed.  The  routine 
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loops  on  the  number  of  body  segments  used  to  specify  the  body 
shape.  Subroutine  B0DPN2  is  called  for  each  of  NFUS  segments  to 
calculate  the  revised  axial  spacing  on  the  body  and  the  body 
panel  geometry.  If  more  than  one  segment  exists,  the  revised 
meridional  spacing  geometry,  YB  and  ZB,  are  read  from  TAPE8.  At 
the  end  all  geometric  arrays  in  blank  common  are  saved  on  TAPE7 . 
The  optional  print  of  source  panel  control  point  coordinates, 
inclination  angles  and  areas  are  performed  here.  A  listing  of 
this  routine  is  found  in  Figure  A-l (e)  of  this  report. 

See  the  discussion  in  Appendix  B,  Section  B-3  regarding 
item  2  in  Figure  A-2  for  the  equivalent  dimension  of  blank  common 
at  this  point.  The  starting  location  indices  defined  in  BODPAN 
and  saved  in  labeled  common  DIMENS  are 

IXPT=1 

IYPT=IXPT+NBODY 
I ZPT=IYPT+NBODY 
ITH=I ZPT+NBODY 
IDEL=ITH+NBODY 
IAR=IDEL+NBODY 

IXC(IFU) =IAR+KFORX ( 1 ) + . . ,+KFORX(IFU)  IFU=1,NFUS 

IYC  ( IFU)  =IAR+KX+KFORX  ( 1)  *KRADX(1)+.  .  .+KFORX(IFU)  'KRADX(IFU) 

IZC (IFU) =IYC ( IFU) +KXKR 

Subroutine  references 

I OREAD ,  B0DPN2,  IOWRIT 

Called  by 
GEOM 

A-8  Subroutine  B0DPN2 

Subroutine  B0DPN2  performs  the  revision  of  the  axial  spacing 
on  noncircular  bodies  and  computes  the  body  panel  georrtetry  of  a 
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single  segment.  The  x,  y,  and  z  coordinates  of  the  cross 
sections  of  a  segment  are  passed  through  the  argument  list  as 
arrays.  A  listing  of  the  routine  is  presented  in  Figure  A-l(f) 
of  this  report. 

The  body  panel  geometry  is  established  by  a  linear  interpola¬ 
tion  along  body  meridian  lines  of  the  y  and  z  coordinates  at  the 
new  axial  stations.  The  interpolation  is  started  with  the  first 
ring  of  panels  in  a  segment  and  continued  until  the  last  ring  of 
panels  on  the  last  segment  is  reached.  The  corner  point  coordi¬ 
nates,  the  control  point  coordinates,  the  inclination  angles,  and 
area  are  calculated  for  each  panel  in  sequence.  The  calculation 
of  the  latter  quantities  are  performed  in  routine  PANEL  from 
corner  point  information. 

The  panel  control  point  coordinates,  the  panel  dihedral  angle 
0,  the  panel  incidence  angle  6,  the  corner  point  coordinates 
and  the  panel  areas  are  returned  through  the  argument  list  to 
BODPAN  where  they  are  stored  in  blank  common.  Finally  if  the 
print  option  IPRT(2)^0,  the  corner  point  coordinates  are 
written  on  the  output  file. 

Argument  list 
JP 

XB  ( M) 

XJ  (J) 

YB (M, K) , 

ZB ( M , K) 

XYZ (JP, 1) , 

XYZ (JP, 2) , 

XYZ (JP, 3) 

XYZ (JP, 4) 


panel  number 

x  station  of  external  geometry  for  Mth  section 
in  segment 

x  station  of  revised  axial  spacing  for  Jth 
section  in  segment 

y  and  z  of  external  geometry  for  Mth  section  in 

segment  and  Kth  meridional  station 

XPT ,  YPT,  and  ZPT  arrays  of  the  x,  y,  and  z 

coordinates  of  the  control  points  of  JPth 

panel 

dihedral  angle  THET  of  the  JPth  panel 
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XYZ(JP,5)  incidence  angle  DELTA  of  the  JPth  panel 

XYZ(JP,6)  area  of  JPth  panel 

XC ( J )  corner  points  of  panels  at  Jth  axial  station 

YC(J,K),  y  and  z  corner  points  of  panels  at  Jth  axial 

ZC ( J ,  K)  station  and  Kth  meridional  angle 

NFUSOR  number  of  external  geometry  stations  on  segment 

KRAD  number  of  meridan  lines  on  segment 

KFUSOR  number  of  panel  revised  axial  stations 

NBODY  total  number  of  body  panels 

IPRT  array  of  print  controls 

Subroutine  references 
PANEL 

Called  by 
BOD PAN 

A-9  Subroutine  BODVEL 

Subroutine  BODVEL  organizes  the  looping  through  source  panels 
to  define  corner  geometry  and  compute  panel  influence  coeffi¬ 
cients.  BODVEL  loops  through  the  number  of  source  panels  on  each 
ring  in  the  given  body  segment  to  set  up  the  local  panel  corner 
coordinates.  A  listing  of  this  routine  is  presented  in  Figure  A-l(f) 
of  this  report. 

The  influence  coefficients  for  each  panel  on  all  remaining 
panels  is  set  up  and  computed  in  PANVEL.  The  influence  coeffi¬ 
cients  are  then  generated  for  one  ring  of  body  panels  on  another 
ring.  For  each  such  block  of  coefficients,  the  influence  normal 
to  the  panel  and  each  of  the  three  component  velocities,  u,  v,  and 
w,  are  saved.  The  block  of  normal  coefficients  is  saved  on  TAPE9. 

If  the  coefficients  correspond  to  a  diagonal  block  in  the  total 
coefficient  matrix  (i.e.  the  influence  of  a  ring  on  itself)  the 
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L*U  tr iangularization  of  the  block  is  performed  by  PAS001  for 
later  calculation  use.  In  supersonic  flow  it  is  assumed  that  aft 
rings  of  panels  have  no  influence  on  upstream  panels.  All  blocks 
of  coefficients  above  the  diagonal  blocks  are  thus  assumed  to  be 
zero  and  not  computed. 

The  descriptions  of  the  variables  in  the  argument  list  follow 


XPT , YPT , ZPT 
THET 
DELTA 
AN 

UB(1, IJ)  , 

UB ( 2 , I J)  , 
UB (3 , IJ) 

XC 

YC  ,  2  C 

KFUSOR 

IFU 

IXZSYM 


coordinates  of  panel  control  points 

array  of  panel  inclination  angles  at  control  points 
array  of  panel  incidence  angles 

temporary  array  to  contain  the  normal  velocities 
to  control  point  I 

arrays  of  velocity  components,  u,v,w,  at 
control  points  in  body  coordinate  system 

x  coordinates  of  leading  and  trailing  edges  of 
panel  rings  of  segment 

y  and  z  coordinates  of  body  panel  corner  points 
of  segment 

number  of  axial  stations  defining  segment 
segment  index 

body  XZ  plane  of  symmetry  indicator 


Subroutine  references 

IOWRIT,  PANVEL,  PAS001 


Called  by 
VELCMP 
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A-l 0  Subroutine  BSHOCK 

Subroutine  BSHOCK  computes  the  nonlinear  shock  wave  shape 
emanating  from  the  store  or  fuselage  nose  produced  by  an 
arbitrary  body  at  angle  of  attack  along  a  given  meridian,  PHIS(J). 
Three  options  are  available  depending  on  the  value  of  NSHOCK 
(input  items  14  or  61)  for  the  choice  of  where  the  meridional 
traverses  are  made.  Single  traverses,  evenly  distributed 
traverses,  and  user  specified  locations  may  be  used  to  define 
polar  shape  of  the  shock  wave.  A  listing  of  the  routine  is 
presented  in  Figure  A-l(g)  of  this  report. 

The  methods  used  to  generate  the  nonlinear  position  of  the 
shock  wave  along  a  single  traverse  are  developed  in  Section  4.1.2 
of  Reference  3.  Use  is  made  of  the  work  of  Reference  2.  The 
axial  scan  for  the  maximum  radial  velocity  behind  the  linear  shock 
wave  is  computed  for  a  fixed  interval  equal  to  half  the  distance 
from  the  nose  to  the  body  shoulder.  This  approximate  interval  must 
be  used  to  account  for  the  saw  tooth  behavior  of  the  velocity  field. 
The  three-dimensional  effects  associated  with  angle  of  attack  are 
accounted  for  by  two  means.  In  computing  the  local  radial 
velocities  along  a  traverse,  the  perturbation  components  must  be 
rotated  into  wind  axes  to  avoid  including  angle  of  attack  components 
directly  in  the  radial  velocities  as  follows 

u  =  u  cosa  -  v  sina  sin^  +  w  sina  cosi|in  +  1.0  (A-2) 

wscsc  Rsc  R 

vw  =  vscos4)r  +  wssin<J)R  (A-3) 

w  =  -u  sina  -  v  cosa  •  sin<t>„  +  w  cosa  •  cos<J>_.  (A-4) 

w  scsc  Rsc  R 

and  second,  the  assumption  the  Mach  cone  rotates  rigidly  about  the 
body  with  angle  of  attack  is  modified  in  accordance  with  refer¬ 
ence  2  as  follows. 
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0 


(A- 5) 


=  6S  +  ac(l  -  ea)  cos<(> 

where  0  is  the  shock  angle  propagating  from  the  nose,  6S  is  the 
nose  limited  shock  angle  at  zero  angle  of  attack,  a  is  the 
included  angle  of  attack,  is  Mach  number  and  nose  cone  angle 
dependent  correction  from  Reference  2,  and  0  is  the  polar  angle 
around  the  cone  of  the  traverse.  Tables  of  RSHK  and  XSHK  for 
each  of  the  polar  traverses  are  the  primary  output  of  this  routine. 

This  routine  is  used  in  two  ways.  For  the  fuselage,  the  shock 
shape  is  computed  at  the  angle  of  attack  of  the  body.  The 
traverses  are  computed  over  the  entire  field  of  interest  and 
include  the  corrections  for  angle  of  attack.  For  the  stores, 
the  shock  shape  is  computed  at  zero  angle  of  attack  for  only  the 
first  quadrant.  The  remaining  shape  is  derived  from  symmetry. 

The  angle  of  attack  correction  of  Equation  (A-5)  is  then  applied 
to  the  zero  angle  of  attack  shape  for  the  local  conditions  of 
each  of  the  elliptic  stores. 

The  descriptions  of  the  parameters  of  the  argument  list  are 
as  follows: 

RINIT  minimum  radius  at  which  shock  shape  search  is 

initiated 

PHISHK  for  NSHOCK=0,  angle  at  which  shock  shape  is 

generated 

for  NSHOCK<0,  maximum  meridional  angle  to  which 
shock  shape  is  calculated  on  even 
spac ing 

THSHK  limiting  shock  wave  angle  at  nose  measured  from 

x-axis,  0g 
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Subroutine  references 
FLDVEL 

Called  by 

WDYBDY,  STORIO 

A-ll  Subroutine  CONFIG 

Subroutine  CONFIG  is  used  to  input  the  geometrical 
description  of  the  external  shape  of  the  noncircular  fuselage  or 
elliptic  stores.  A  listing  of  the  routine  is  presented  in  Figure 
A-l(i)  of  this  report.  The  routine  first  reads  the  configuration 
reference  area  from  the  input  file  if  J0/0,  otherwise  the 
reference  area  is  set  equal  to  unity. 

If  J2/0,  the  body  external  geometry  data  is  read  from  the  input 
file  according  to  the  specified  body  option  for  each  of  NFUS  body 
segments.  For  arbitrary  cross-sections,  the  y  and  z  ordinates  of  the 
body  segment  are  read.  For  circular  or  elliptic  sections,  the 
particular  combination  of  body  area,  radius,  semi-major  or  minor 
axes,  or  elliptic  ratio  requested  may  be  read. 

In  addition,  the  area  at  each  axial  station  is  computed.  The 
maximum  cross-section  area  is  also  saved.  The  description  of  the 
single  argument  is  as  follows: 

SFUS  array  to  hold  the  y  and  z  ordinates  of  the  arbitrary 

cross-section  data  for  J2=l 

Subroutine  references 
ASECTN 

Called  by 
GEOM 
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A-12  Subroutine  DOUBLT 


Subroutine  DOUBLT  calculates  coefficients  used  in  the 
determination  of  the  line  doublet  strengths.  They  occur  as 
terms  in  Equations  (1-27)  and  (1-28)  of  Reference  1.  The 
relation  of  the  coefficients  to  perturbation  velocities, 
uB  and  vB  ^/V^,  induced  by  a  number  of  line  doublets 

distributed  along  the  body  centerline  is  shown  in  the  first 
two  of  Equation  (1-30)  in  which  the  coefficients  occur  as  the 
multipliers  of  Kd  ncos  0.  A  listing  of  the  subroutine  is 
presented  in  Figure  A-l(j)  of  this  report.  Section  1-2.2  of 
Reference  1  should  be  referred  to  for  further  details  concerning 
the  doublet  strength  calcula t ions . 

For  a  specified  control  point,  x  ,  rD,  and  singularity  origin 

D  D 

,  the  subroutine  calculates  quantities  U  and  V  according  to  the 
following  equations. 


V 


cosh 


rxB  " 

1  +  ; 

fXB  '  M 

\ 

/ 

fXB  ' 

>  ®rB  • 

1  + 

®rB  • 

l  6rB  J 

( A-6 ) 


(A-7) 


At  the  beginning  of  the  subroutine  a  test  is  performed  to  deter¬ 
mine  if  the  control  point  is  ahead  of  the  Mach  cone  from  the 
doublet  origin.  If  so,  U  and  V  are  set  to  zero  and  control  is 
returned  to  the  calling  program. 

The  following  table  of  definitions  contains  most  of  the 
variable  names  used  in  the  subroutine: 
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BUT  A 

RFIELO 


XF1ELD 


r^;  radius  of  body  at  control  point 

;  location  on  body  axis  of  doublet  origin, 

positive  measured  from  tip  of  nose 

coefficient  defined  by  Equation  A-6 

coefficient  defined  by  Equation  A-7 

x  ;  x  location  of  control  point,  positive 
B 

measured  from  tip  of  nose 


Called  by 
BDYGEN 


A-13  Subroutine  DPCOEF 


Subroutine  DPCOEF  calculates  the  coefficient  matrix  of  the 
set  of  simultaneous  boundary  condition  equations  which  are  to  be 
solved  for  the  constant  u-velocity  panel  singularity  strengths. 

A  listing  of  the  subroutine  is  presented  in  Figure  A-l(k),  and  a  flow 
chart  in  Figure  A-5  of  this  report. 

The  elements  of  the  matrix  are  the  aerodynamic  influence 
coefficients  described  in  Section  3.3.4  of  Reference  3.  They 
occur  indirectly  in  the  summation  terms  in  the  left-hand  side 
of  Equations  (24),  (25),  and  (26)  of  that  reference.  The  actual 

coefficients  which  the  subroutine  calculates,  FVN(v,n),  are 
related  to  the  summation  terms  through  the  panel  strengths, 
u^/V^,  by  the  following  equations  in  which  n  is  the  index  of 
the  influencing  panel  and  v  is  the  control  point  index: 


cos4> 


FVN ( v , n ) 


N  u 

v, n  _  +n 


FVN ( v , n ) 


FVN (v,n) 


v  =  1,2,.. . NPANLS 


v  =  NIP ,NIP+1 , . . . N2 


v  =  N2P,N2P+1, . . .NPTOT 


for  n  =  1,2,... , NPTOT . 

The  subroutine  consists  of  three  double  DO  loops.  The  first 
loop  uses  subroutine  VELWP1  to  calculate  the  influence  of  the  wing 
panels  at  the  wing,  pylon,  and  fuselage  control  points.  The  second 
loop  is  bypassed  if  there  is  no  pylon,  NPY=0.  If  e  pylon  is 
present,  the  influence  coefficients  for  the  pylon  panels  are 
calculated  by  means  of  subroutine  VELPP1.  The  third  double  loop 
is  bypassed  if  there  is  no  fuselage,  NFU=0.  Otherwise,  the 
influence  coefficients  for  the  fuselage  panels  are  calculated 
using  subroutine  VELBD1. 

Within  each  outer  loop  are  three  inner  loops  in  series,  which 
fix  the  control  point  location  on  the  wing,  pylon,  or  fuselage, 
respectively.  On  the  flow  chart,  only  the  first  occurrence  of 
the  first  inner  loop  is  shown  in  detail,  which  includes  the  call 
to  subroutine  VELWP1.  The  remaining  inner  loops  have  the  same 
logical  structure  and  are  not  shown. 

If  the  control  point  is  located  on  a  wing  or  a  fuselage  panel, 
the  influence  coefficient  is  the  component  normal  to  the  panel 
surface  at  the  control  point.  Thus,  the  wing  coefficients  are 
rotated  through  the  panel  dihedral  angle  using  quantities  pre¬ 
viously  calculated  in  subroutine  WLYOUT.  The  fuselage  coeffi¬ 
cients  are  rotated  through  the  panel  orientation  angle  using 
quantities  calculated  in  subroutine  BLYOUT. 

Subroutine  references 

VELBD1,  VELPP1,  VELWP1 

Called  by 

LDCALC 


A-14  Subroutine  DPRHS 


Subroutine  DPRHS  calculates  the  right-hand  side  vector  of  the 
set  of  simultaneous  boundary  condition  equations  which  are  to  be 
solved  in  order  to  determine  the  constant  u-velocity  panel 
singularity  strengths.  The  boundary  conditions  are  specified 
in  Equations  (24),  (25),  and  (26)  of  Reference  3.  A  listing  of 
the  subroutine  is  presented  in  Figure  A-l(k)  and  a  flow  chart  in 
Figure  A-6  of  this  report. 

The  major  portion  of  the  subroutine  is  devoted  to  evaluating 
the  externally  induced  perturbation  velocities,  uw  /V  , 

w  1  '  v  OO 

V  /Vm ,  and  ww .  /VM,  at  all  of  the  wing,  pylon,  and  fuselage 

1  §  V  1  (  V 

control  points. 

The  first  section  calculates  velocities  induced  at  wing  and 
pylon  control  points  by  the  circular  fuselage  line  sources,  sinks, 
and  line  doublets  if  a  circular  fuselage  is  present  (NFU=1).  A 
control  point  is  located  in  the  fuselage  coordinate  system  and 
subroutine  VEIXAL  is  called  to  calculate  the  velocities  at  this 
poin*-.  These  velocities  are  summed  in  the  UEI,  VEI,  and  WEI 
a  r rays . 

The  next  section  calculates  velocities  induced  at  wing  and 
pylon  control  points  by  the  noncircular  fuselage  source  panels  if 
a  noncircular  fuselage  is  present  (NFU=2)  .  Data  for  the  fuselage 
are  first  restored  to  blank  common  from  TAPE10  through  FRSTRT  and 
indices  defining  the  beginning  of  the  various  arrays  are  defined. 
All  of  the  wing  and  pylon  control  points  are  then  located  in  the 
fuselage  coordinate  system.  The  velocities  at  these  points  are 
calculated  by  calling  subroutine  FLDVEL. 

Next,  the  perturbation  velocities  induced  by  the  wing  and 
pylon  thickness  source  panels  are  calculated  by  means  of 
subroutines  VELWT1  and  VELPT1.  The  effects  included  are  wing 


source  panels  on  wing,  pylon;  and  fuselage  control  points  and 
pylon  source  panels  on  wing  and  fuselage  control  points.  If  a 
pylon  is  not  present,  NPY=0,  or  the  fuselage  is  not  present, 

NFU  =  0 ,  calculations  involving  the  nylon  or  the  fuselage  are 
omitted. 

The  last  three  loops  in  the  subroutine  calculates  the 
right-hand  sides  of  the  equations.  If  the  pylon  is  located  below 
the  fuselage  centerline,  CENTER=TRUE,  VE I  is  set  to  zero  for  all 
pylon  control  points. 

It  should  be  noted  that  if  the  control  point  is  located  on  the 
wing  or  on  the  f  iselage,  dihedral  effects  are  included  in  calcula¬ 
ting  the  velocity  normal  to  the  panel  surface. 

Subroutine  references 

FLDVEL,  FRSTRT,  VELCAL,  VEIPT1,  VELWT1 

Called  by 
LDC  ALC 

A- 15  Subroutine  FLDVEL 

Subroutine  FLDVEL  is  used  to  organize  the  computation  of  the 
u,  v,  w  velocity  components  at  the  field  points  XFP,  YFP,  ZFP. 

This  routine  initializes  all  velocities  to  zero  and  computes  the 
locations  of  the  required  geometric  and  strength  arrays  in  blank 
common  for  the  fuselage  or  store  at  hand.  If  the  configuration 
models  the  fuselage  with  inlets,  variable  INLET  is  initialized  to 
TRUE  to  indicate  presence  of  inlet  panels.  It  then  performs  the 
looping  for  the  number  of  segments  and  the  number  of  rings  in  each 
segment  to  sum  the  component  velocity  contributions  at  each  field 
point.  A  listing  of  the  routine  is  presented  in  Figure  A-l(m)  of 
this  report. 


The  descriptions  of  the  parameters  in  the  argument  list 

to  1  low : 

XFP,YFP,ZFP  arrays  of  coordinates  of  field  points  in  coordinate 
system  of  body  panels  at  which  component  velocities 
are  computed 

SVN  temporary  array  of  length  NFLD  used  to  test  for 

zero  influence  at  field  point 

LSKP  temporary  logical  array  of  length  NFLD  used  to  test 

for  zero  influence  of  ring  of  panels  at  field  point 

U,V,W  arrays  of  component  velocities  computed  for 

influence  of  body  at  field  points 

NFLD  number  of  field  points 

ID  array  containing  the  same  variables  in  the  same  order 

as  in  common  DIMENS 

IG  index  specifying  IGth  strength  solution,  GB,  to 

be  used 

Subroutine  references 
FLDVL2 

la  1  led  by 

BSHOCK,  DPRHS,  INLSHK 

A-16  Subroutine  FLDVL 2 


Subroutine  FLDVL2  computes  the  three  components  of  velocity 
induced  at  specified  field  points  by  the  source  panels  of  a  given 
ring  of  body  panels.  The  field  point  is  assumed  to  have  no 
incidence  or  inclination  relative  to  the  flow.  A  listing  of  the 
routine  is  presented  in  Figure  A-l  (m)  of  this  report. 

The  routine  first  initializes  the  field  point  inclination  and 
incidence  transformations  to  zero.  For  each  of  the  influencing 


panels  on  a  ring  the  control  point  coordinates,  inclination  and 
incidence  angles  and  corner  points  in  the  local  panel  system  are 
defined.  For  each  panel  g^  is  initialized  to  g  .  For  panels  used 
to  model  inlet  openings  g^  is  set  to  the  minimum  of  g  and  BTINLT. 
Further,  if  an  inlet  panel  slope  exceeds  BTINLT,  g^  is  set  to 
0.99/tan6.  For  each  of  the  field  points  the  influence  of  the 
panel  on  the  field  point  is  computed,  multiplied  times  the  panel 
strength  and  summed. 

To  minimize  computational  time  FLDVL2  tests  whether  each  ring 
of  panels  contributes  to  the  total  velocity  at  a  given  field  point. 
For  each  point  the  following  sum  is  made  for  each  panel  on  the  ring 

SVN  (I)  =  UB2  +  VB2  +  WB2 


where  JB,  VB ,  WB  are  the  influence  coefficients  of  a  panel  at  the 
Ith  field  point.  If  the  net  influence  for  a  complete  ring  SVN  is 
equal  to  zero  the  logical  variable  LSKP(I)  for  that  field  point  is 
set  .TRUE,  and  all  further  calculations  at  that  point  are  sup¬ 
pressed.  This  is  based  on  the  assumption  that  at  supersonic  speeds 
once  a  point  is  ahead  of  the  Mach  waves  from  a  ring,  the  influence 
of  that  ring  and  all  subsequent  rings  will  vanish. 

The  descriptions  of  the  parameters  in  the  argument  list  follow: 


XPT, YPT, ZPT 

THET 

DELTA 

GB 

XFP, YFP, ZFP 

u,v,w 


arrays  of  the  coordinates  of  the  source  panel 
control  points  for  the  entire  body 
array  of  inclination  angles  for  panels 
array  of  incidence  angles  for  panels 
array  of  panel  strengths 
arrays  of  coordinates  of  field  points 
arrays  of  orthogonal  velocity  components  at 
panels  in  x,  y,  and  z  directions 
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SVN 


temporary  array  used  to  test  influence  of  ring 
on  po  i n  t 

LSKP  temporary  array  used  as  logical  indicator  of  no 

further  influence  at  point 
NFLD  number  of  field  points 

XC , YC , ZC  arrays  of  coordinates  of  panel  corners  for 

segment 

KFUSOR  number  of  axial  station  bounding  panels  in 

segment 

IXZSYM  body  symmetry  indicator 

JR  number  of  panels  in  ring  J 

JG  starting  offset  index  of  first  panel  in  ring 

L  index  of  trailing  x-station  in  segment 

Subroutine  references 
PANVEL,  INLTST 

Called  by 

FLDVEL,  IMAGEV 

A-l 7  Subroutine  FRSTRT 

Subroutine  FRSTRT  is  used  to  save  or  restore  required  program 
information  for  the  source  paneling  method  to  restart  a  con¬ 
figuration  analysis.  The  routine  is  set  up  to  allow  storage  of 
the  information  from  more  than  one  configuration  on  the  same 
file.  A  listing  of  the  routine  is  presented  in  Figure  A-l{n)  of 
this  report.  Six  options  are  available  based  on  the  information 
to  be  stored  or  retrieved  and  where  the  information  is  to  reside. 
They  and  the  functions  they  perform  are 

DESCRIPTION 

saves  control  integers  and  arrays  in  blank  common 
restores  control  integers  and  arrays  in  blank  common 


KODE 

1 

2 


3 


saves  above  information  and  AIC,  U,  V,  and  W 
velocity  matrices 

4  restores  above  information  and  AIC,  U,  V,  and  W 

velocity  matrices 

5  restores  control  integers  and  arrays  stored 

under  either  option  1  or  3  into  starting  core 
location  in  blank  common.  All  arrays  are  copied 
into  blank  common  rather  than  into  original 
labeled  commons 

6  restores  information  in  manner  of  K0DE=5  and  reads 

past  velocity  matrices  stored  under  option  3  to 
position  file  at  next  record 

The  information  saved  under  K0DE=1  consists  of  all  the  varia¬ 
bles  in  labeled  commons  DIMENS,  PARAM,  BOPTNS,  BGEOM,  HEAD,  and 
BSHOCK,  and  the  first  NTAP7  variables  in  blank  common.  The  later 
set  are  also  saved  on  TAPE7 .  If  the  configuration  contains  an 
inlet,  the  variables  in  common  BINLET  and  BINSHK  are  also 
saved.  In  copying  the  velocity  matrices  under  K0DE=2,  the  arrays 
are  first  copied  from  TAPE8  or  TAPE9  into  temporary  arrays  in 
blank  common  and  then  onto  the  output  file  10.  The  last  two 
options  are  used  in  Program  II  to  stack  multiple  configurations 
end-to-end  in  blank  common.  Provision  has  been  made  for  only  one 
configuration  with  an  inlet.  Variables  from  commons  BINLET  and 
BINSHK  are  copied  back  into  those  arrays  under  all  KODE  options. 
During  this  data  retrieval  three  indices  used  to  locate  the 
control  variable  arrays  are  saved. 

IDO  (=LASTA+1)  first  location  in  blank  common  of  variables 

contained  in  labeled  commons  DIMENS,  PARAM, 
BOPTNS,  BGEOM,  and  TITLE  in  order 

ISKO  (=IDO+571)  first  location  in  blank  common  of  variables 

contained  in  labeled  common  BSHOCK 
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I AO  (=ISKO+240) 


first  location  in  blank  common  of  geometric  and 
strength  arrays  previously  saved  on  TAPE7 


The  descriptions  of  the  parameters  in  the  argument  list  are: 

10  external  file  unit  number  onto  which  the  data 

is  stored  or  retrieved  in  unformatted  form 

KODE  optional  index  selecting  information  to  be  read 

or  written  to  file  10;  see  above  descriptions 

LASTA  last  location  in  blank  common  currently  defined 

Subroutine  references 
I OREAD ,  IOWRIT 

Called  by 

WDYBDY,  STORIO,  DPRHS ,  WRFILE 

A-18  Subroutine  FUSEIO 

Subroutine  FUSEIO  reads  and  prints  the  input  data  which 
describe  the  circular  fuselage  (NFU=1)  and  calculates  the  line 
source  and  doublet  distributions  as  described  in  Appendix  I  of 
Reference  1.  A  listing  of  the  subroutine  is  presented  in 
Figure  A-l(o)  and  a  flow  chart  in  Figure  A- 7  of  this  report. 

The  subroutine  first  reads  in  and  prints  input  items  5,  6,  7, 
and  8  which  consist  of  the  fuselage  length,  the  maximum  radius, 
and  the  polynomials  specifying  the  fuselage  shape.  Next,  the  data 
used  to  lay  out  the  body  interference  panels,  input  items  9  and 
10,  are  read  and  printed.  Subroutine  BDYGEN  is  then  called  for 
the  purpose  of  calculating  the  source  and  doublet  distributions. 
The  fuselage  shoulder  location  is  next  found  and  finally 
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subroutine  SHKSHP  is  called  to  calculate  the  nose  shock  wave 
shape . 

Subroutine  references 
BDYGEN,  SHKSHP 

Called  by 
LDCALC 

A-19  Subroutine  GEOM 

Subroutine  GEOM  organizes  the  reading  of  the  noncircular  body 
input  and  the  calculation  of  the  source  panel  geometry  arrays. 

This  routine  is  used  to  input  both  the  noncircular  fuselage  and 
the  elliptic  store  panel  geometry.  A  listing  of  the  routine  is 
presented  in  Figure  A-l(p)  of  this  report.  * 

The  specification  of  the  input  describing  the  body  paneling 
is  split  in  two  phases  by  GEOM.  In  the  first  phase  the  control 
parameters  and  the  geometric  arrays  describing  the  external  shape 
of  the  body  are  read.  The  external  geometry  arrays  are  read  by  a 
call  to  CONFIG.  This  geometry  defines  the  external  shape 
independently  of  the  paneling  used.  It  may  optionally  be  used  to 
specify  panel  corner  coordinates.  In  the  second  phase  the  control 
parameters,  additional  arrays  and  inlet  data  required  are  read 
which  define  the  subdivision  of  the  external  shape  into  panel 
coordinates.  The  routine  NEWRAD  is  called  to  perform  the 
interpolation  necessary  to  redefine  the  meridional  spacing  of  the 
panels.  The  routine  BODPAN  is  called  to  perform  the  inter¬ 
polations  for  new  axial  spacings  of  panels  and  to  calculate 
control  point  coordinate  and  inclination  information. 

This  routine  also  computes  the  number  of  panels  to  be  generated 
and  the  allocation  of  array  space  in  blank  common  to  be  used. 

The  total  number  of  axial  stations,  meridional  angles  used  to 


define  panel  corners  and  the  total  length  of  blank  common 
required  are  compared  against  dimensioned  values  as  check  on 
input  problem  size. 

Subroutine  references 

BODPAN,  CONFIG,  NEWRAD 

Called  by 

WDYBDY,  STORIO 

A-20  Subroutine  INLSHK 

Subroutine  INLSHK  computes  the  nonlinear  shock  wave  shape 
emanating  from  a  ramp  type  inlet  shown  in  Figure  12  of  Volume  II. 
The  table  of  shock  locations  is  produced  in  a  manner  analygous  to 
that  used  in  BSHOCK  to  generate  the  shock  shape  about  the  body 
nose  with  the  exceptions  that  only  the  influence  of  the  modified 
panels  across  the  inlet  are  felt  and  the  shock  propagates  from  the 
leading  edge  of  the  ramp.  The  detailed  explanation  of  the  inlet 
shock  is  presented  in  Section  4.4.2  of  Reference  3.  A  listing  of 
the  routine  is  presented  in  Figure  A-l(r)  of  this  report. 

The  generation  of  the  traverses  used  to  locate  the  first 
influence  and  integrate  the  nonlinear  shape  uses  the  same  model 
and  equations  detailed  for  BSHOCK.  Additional  assumptions  are 
made  about  the  body  in  generating  the  traverses  regarding  the 
shape  of  the  ramp  inlet.  Field  points  aft  of  the  inlet  are 
assumed  to  lie  in  the  three  regions  relative  to  XINLT,  YINLT, 

ZINLT  as  shown  in  the  following  sketch. 


REGION  II 

pIfs(y’z) 


In  region  I  inside  and  below  the  inlet,  traverses  are  carried  out 
parallel  to  the  XZ  plane  with  subsequent  interpolation  in  the 
y-direction  to  find  the  first  inlet  shock  influence,  g  .  in 

.  1  I  S 

region  II,  outboard  of  the  inlet,  traverses  are  calculated  along 
radial  lines  from  (XINLT,  YINLT,  ZINLT)  from  4)  =  0°  to  <j>  =  <J> 

SHK 

The  last  traverse  at  4>  =  180°  is  assumed  such  that  8  =  8  .  In 

I ,  S  00 

region  III  inboard  and  above  4>g  for  the  inlet,  the  flow  is  assumed 

equal  to  free  stream  for  which  6  is  set  to  p. 

I ,  S 

NIS  shock  traverses  are  generated  alternating  first  in  Region  I 
and  then  in  Region  II  to  define  the  inlet  shock  location.  The 
shocks  propagate  from  a  starting  point  at  (XINLT,  y,  ZINLT)  along 
a  line  at  the  angle  1  from  the  vertical  plane.  The  first  shock 
is  positioned  at  y=YINLT  and  4'  =  0°.  The  second  traverse  propagates 


from  y=YINLT  at  <)  =  PHISHK.  The  third  propagates  from  y=0  at  <J>  =  0°. 
The  fourth  is  at  y=YINLT  and  ,^=4  5°  and  the  fifth  is  at  y=YCPI  and 
;=0° .  Additional  traverses  are  laid  out  by  subdividing  region  I 
between  y=0  and  y=YCPI  and  region  II  between  i>  =  0  and  <j)=PHISHK  with 
even  increments. 

The  following  are  differences  in  the  calculative  procedure  for 
the  inlet  and  the  nose  shock  shape.  In  scanning  in  the  axial 
direction  at  a  given  Y,Z  location,  all  values  are  referenced  to 
XINLT ,  YINLT(IRAD),  and  ZINLT.  In  computing  the  influence  at 
field  points  the  strengths  of  all  panels  except  those  designated 
as  open  or  blocked  inlet  panels  and  additional  adjacent  panels 
on  the  inlet  aft  of  the  opening  are  set  equal  to  zero.  In 
integrating  the  traverses,  no  angle  equivalent  to  the  nose 
limiting  shock  angle,  6  ,  is  used.  When  estimating  the  shock 
location  below  the  inlet  in  region  I  in  the  presence  of  partial 
blockage  (mass  flow  ratios  less  than  one) ,  the  shock  is  assumed 
to  follow  BTINLT  down  the  inlet  face  to  the  first  blocked  panel. 

A  normal  shock  is  assumed  to  exist  moving  the  shock  vertically 
down  to  the  lower  lip  of  the  cowl.  For  the  first  step,  thereafter, 
in  the  radial  direction,  the  6  c  propagating  from  the  inlet  is 
set  equal  to  the  first  computed  value.  The  only  restriction  is 
that  it  be  positive. 

The  description  of  the  parameter  in  the  argument  list 
follows : 

PHISHK  maximum  meridional  angle  about  YINLT, ZINLT  to 

which  shock  shape  is  calculated  on  even  spacing 

Subroutine  references 

FLDVCL,  INLTST 


Called  by 
WDYBDY 
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A-21  Function  INLTST 

Logical  function  INLTST  determines  whether  the  panel  index  in 
the  argument  is  for  an  inlet  panel.  The  value  of  the  function  is 
set  to  TRUE  if  the  index  is  for  an  inlet  panel  and  FALSE  if  it  is 
not.  In  addition  the  logical  variable  OPEN  is  set  to  indicate 
whether  the  inlet  panel  is  blocked  or  unblocked  to  flow.  Three 
conditions  are  tested  for.  If  no  inlet  panels  exist  INLTST  and 
OPEN  are  set  false.  If  an  inlet  panel  exists  the  index  I  is 
compared  to  the  table  of  possible  inlet  panel  numbers,  JINLT.  If 
I  is  equal  to  one  of  the  inlet  panel  numbers,  INLTST  is  set  true. 

If  I  is  not  an  inlet  panel  number, INLTST  is  set  false.  If  I  falls 
in  the  subset  of  JINLT  of  unblocked  panels,  OPEN  is  true; 
otherwise  OPEN  is  set  false.  A  listing  of  this  routine  is 
presented  in  Figure  A-l(t)  of  this  report. 

The  description  of  the  parameters  in  the  argument  list  follows. 

I  panel  number  index  to  be  compared  with  table  of 

possible  inlet  panel  numbers,  JINLT 
OPEN  logical  variable  indicating  whether  an  inlet  panel 

is  open  or  blocked.  OPEN  is  TRUE  if  panel  allows 
unblocked  flow  through  panel 

Called  by 

BODVEL,  FLDVL2 ,  PANEL,  SOLVE 


A-22  Subroutine  INLXYZ 

Subroutine  INLXYZ  scans  the  X,Y,Z  coordinates  of  the  inlet 
panel  corners  to  find  the  coordinates  to  be  used  to  define  inlet 
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leading  edge,  lower  cowl  lip,  and  first  axial  station  containing 

blocked  panels.  Nine  values  are  computed.  They  correspond  to  the 

most  outboard  leading  edge  of  the  inlet  panels,  XINLT, YINLT, ZINLT, 

the  most  aft  outboard  trailing  edge  of  the  inlet  panels,  XINLTE, 

YINLTE, ZINLTE,  the  innermost  y  station  of  the  inlet,  YCPI,  and  the 

first  axial  station  of  blocked  inlet  panels,  XCLOSD.  If  no  inlet 

blockage  exists, this  value  is  set  to  inlet  trailing  edge  XINLTE. 

In  addition,  the  inlet  mass  flow  ratio  is  computed  from  the  ratio 

of  exposed  open  frontal  area  of  inlet  panels,  A  ,  to  total 

inlet  frontal  area,  Atota^«  The  mass  flow  ratio  RVIVO  is  set 

equal  to  A  „  /A.  .  ,.  This  routine  is  called  once  for  each  inlet 

open  total 

panel.  A  listing  of  the  routine  is  presented  in  Figure  A-l(t)  of 
this  report. 


The  descriptions  of  the  parameters  in  the  argument  list  follow. 


XIN, YIN, ZIN 

FAREA 

OPEN 


coordinates  of  corner  of  inlet  panel  to  be 
checked 

frontal  area  of  panel 

logical  variable  indicating  whether  an  inlet  panel 
is  open  or  blocked.  OPEN  is  TRUE  if  panel  allows 
unblocked  flow  through  panel 


Called  by 
PANEL 


A-23  Subroutine  INVER1 

Subroutine  INVER1  solves  the  system  of  simultaneous  linear 
algebraic  equations. 

[A] X  =  B 

This  routine  performs  pivot  searching  during  the  solution  of  the 
general  matrix,  A.  The  right  hand  side  is  passed  int'-*  INVER1  as 
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the  N+l'st  column  in  the  matrix.  The  solution,  X,  also  returns 
in  that  location.  A  listing  of  this  routine  is  presented  in 
Figure  A-l(u)  of  this  report.  The  routine  is  currently  limited  to 
200  equations  by  internal  dimensions.  The  descriptions  of  the 
parameters  in  the  argument  list  follow. 


A 


NS  YS 
N 

NMAX 

MMAX 


coefficients  of  linear  system  of  equations  in  first 

N  columns;  columns  N+l  through  N+NSYS  contain 

multiple  right-hand  sides,  B,  on  input  and 

solutions,  X,  on  return 

number  of  right  hand  sides 

actual  number  of  equations 

first  dimension  of  A 

second  dimension  of  A 


Called  by 
LDCALC 


A-24  Subroutine  IOREAD 

Subroutine  IOREAD  performs  an  unformatted  read  from  external 
file,  10.  NA  consecutive  elements  of  array,  A,  are  read 
sequentially.  This  routine  is  used  to  specify  a  common 
interface  to  external  files.  A  listing  of  tins  routine  is 
presented  in  Figure  A-l(u)  of  this  report.  The  descriptions  of 
the  parameters  in  the  argument  list  follow. 

10  external  file  reference  number 

A  array  of  numbers  to  be  read 

NA  number  of  elements  of  A  to  be  read 

A  machine  dependent  version  of  IOREAD  is  available  for  CDC 
machines  using  input  routine  BUFFER  IN  by  appropriate  replace¬ 
ment  of  comment  cards  within  the  routine. 
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Called  by 

BODPAN,  VELCMP,  SOLVE,  SMARCH,  FRSTRT 

A-25  Subroutine  IOWRIT 

Subroutine  IOWRIT  performs  an  unformatted  write  to  external 
file,  10.  NA  consecutive  elements  of  array,  A,  are  written 
sequentially.  This  routine  is  used  as  a  common  interface  to 
external  files.  A  listing  of  this  routine  is  presented  in 
Figure  A-l(u)  of  this  report.  The  descriptions  of  the  parameters 
in  the  argument  list  follow: 

10  external  file  reference  number 

A  array  of  numbers  to  be  written 

NA  number  of  elements  of  A  to  be  written 

A  machine  dependent  version  of  IOWRIT  is  available  for  CDC 
machines  using  output  routine  BUFFER  OUT  by  appropriate 
replacement  of  comment  cards  within  the  routine. 

Called  by 

BODPAN,  VELCMP,  SOLVE,  NEWRAD,  FRSTRT,  BODVEL 

A-26  Subroutine  NETCLC 

Subroutine  NETCLC  calculates  the  net  corner  strengths  for  a 
given  set  of  wing,  pylon,  or  fuselage  constant  u-velocity  or 
thickness  panels.  The  net  strength  for  a  corner  point  is 
calculated  by  superposing  the  strengths  of  the  four  panels 
surrounding  the  point.  Whenever  a  corner  point  lies  on  an 
exterior  boundary  or  along  a  chordwise  row  where  a  break  in  sweep 
angle  or  dihedral  angle  occurs,  surrounding  panels  are  assumed  to 
be  present  with  zero  strengths.  The  superposition  scheme  is 
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described  in  Section  3.3.7  of  Reference  3.  A  listing  of  the  program 
is  presented  in  Figure  A-i tv)  of  this  report. 

The  panel  strengths  are  passed  to  the  subroutine  as  formal 
parameters  in  the  array  S,  which  is  dimensioned  NCW  by  MSW .  The 
subroutine  first  stores  this  array  in  a  new  array  P,  leaving  the 
first  and  last  columns  zero  and  inserting  a  zero  column  whenever 
a  break  in  leading-edge  or  t r a i 1 i ng -edg e  sweep  angle  or  dihedral 
angle  occurs.  The  polar  angle  acts  as  dihedral  angle  for  a  body 
interference  panel.  The  first  and  last  rows  of  P  are  also  set  to 
zero.  These  zero  elements  of  array  P  correspond  to  panels  of  zero 
strength  mentioned  above. 

The  rest  of  the  subroutine  consists  of  a  double  DO  loop  in 
which  the  net  corner  strengths,  DP,  are  calculated  according  to 
the  following  relation: 

DP  (I,  J)=P  (I,  J)-P  (1  +  1,  J)-P  (I,  J+D+P  (1  +  1,  J  +  l)  ,  ( A—  8 ) 

1=1,  NCW+1 
J=l,  MSW+k+1 

where  k  is  the  number  of  breaks  in  sweep  or  dihedral  which  occur 
over  the  specified  set  of  panels. 

The  descriptions  of  the  parameters  in  the  argument  list  follow: 

S  array  containing  the  NCW*MSW  panel  strengths 

DP  array  containing  the  NCW1*MSW1  net  corner  strengths 

NCW  number  of  panels  in  a  chordwise  row;  also  row 

dimension  of  S 

NCW1  NCW+1;  row  dimension  of  DP 

MSW  number  of  panels  in  a  spanwise  row  on  wing  or 

pylon;  number  of  panels  in  a  ring  of  fuselage  panels 
MSW 1  MSW+1 
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m 


MSWK  1 


PS  ILE 

PS  ITE 


PHI 


number  ot  corner  points  in  a  s panwise  row  of  wing 
or  pylcn  panels;  number  of  corner  points  in  a  ring 
of  fuselage  panels 

array  containing  leading-edge  sweep  angles  for 
wing  or  pylon,  polar  angles  for  fuselage  panels 
array  containing  t r a i 1 ing -edg e  sweep  angles 
for  wing  or  pylon,  polar  angles  for  fuselage 
panels 

array  containing  dihedral  angles  for  wing,  polar 
angles  for  fuselage;  zero  array  for  pylon 


Called  by 
NULYT 


A-27  Subroutine  NEWRAD 

Subroutine  NEWRAD  is  used  to  organize  the  revision  of  the 
noncircular  body  meridian  line  spacing  for  the  source  panels. 

This  routine  computes  the  number  of  meridional  lines  to  be  used 
for  paneling  on  the  half  or  full  body.  The  allocation  of  arrays 
to  be  used  for  temporary  storage,  YB  and  ZB,  are  defined  also. 
NEWRAD  calls  NEWRD2  for  the  revised  meridional  spacings  for  each 
of  the  NFUS  body  segments.  If  more  than  one  segment  exists,  the 
revised  ordinates,  YB  and  ZB,  are  written  on  external  file 
number  8.  A  listing  of  this  routine  is  presented  in  Figure  A-l(v) 
of  this  report. 

Subroutine  references 
IOWRIT,  NEW RD2 


A-28  Subroutine  NEWRD2 

Subroutine  NEWRD2  computes  the  revised  noncircular  body 
meridional  line  spacing  of  a  given  body  segment.  For  each 
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segment,  there  are  three  options  for  redefining  t  tie  meridian 

lines  and  the  option  to  define  the  symmetric  portion  of  ^  nod y. 

It  KRAD=0,  the  meridian  lines  are  not  changed.  If  KRAD  is 

positive,  the  meridian  lines  are  relocated  at  KRAD  equally 

spaced  values  of  the  meridian  angle  •  .  If  KRAD  is  negative, 

the  user  defined  values  of  are  used.  Fo  i  each  of  the  above 

k 

options,  new  meridian  lines  may  be  defined  from  circular, 
elliptic,  or  arbitrary  section  input.  A  listing  of  this 
routine  is  presented  in  Figure  A-l(v)  of  this  report. 

If  the  body  has  circular  or  elliptic  cross  section,  the  y 
and  z  coordinates  are  calculated  at  each  station  as  follows: 


z  =  zc  +  r  sin  : 

where  z c  is  the  camber  offset  in  item  19  of  the  input  and  trie 
local  radius  r  is  defined  as 

r  ~  .  — . — — - - -  i  A  -  i 

/("cos b /BY )  2  +  (sin;/AZ)  2 

The  parameters  BY  and  AZ  are  the  horizontal  and  vertical  semi¬ 
axis  of  the  ellipse.  For  the  circular  body,  BY  and  AZ  are 
equal  to  the  radius  and  the  above  equation  reduces  to  the  circular 
radius. 


If  the  body  has  an  arbitrary  cross  section,  the  y  and  z  coordi¬ 
nates  are  obtained  by  linear  interpolation  at  the  new  values  of 
of  the  original  y  and  z  coordinates  read  in  the  program  input 
i terns  20  and  21. 

If  the  symmetry  option  is  used,  IXZSYM=1,  the  remaining 
coordinates  for  the  second  half  of  the  body  are  defined.  The 
meridian  lines  for  the  symmetric  half  of  the  body  continue  in 
counterclockwise  fashion  about  the  body.  The  y  ordinates  are  t he 
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the  first  half  of  the  body,  while  the  z  values  remain 


iv'oativo  of 
ti.i'  same. 

The  descriptions  of  the  parameters  in  the  argument  list  follow: 

YB,ZB  y  and  z  ordinates  of  the  revised  meridional 

spac  i  ng 

PHIK(K)  ,  Kth  revised  meridian  for  new  spacing,  deg. 

ZFUS(N)  z  ,  camber  offset  of  centerline  at  Nth  axial 

station  of  segment 

FUSBY(N),  BY  and  AZ ,  horizontal  and  vertical  semi-axes  of 

FUSAZ(N)  ellipse  at  Nth  axial  station  of  segment 

bFUS  array  containing  y  and  z  for  arbitrary  cross 

section  input  option 

NFUSOR  number  of  axial  stations  of  segment 

KRAD  number  of  meridian  lines  of  revised  spacing 

NRAD  number  of  meridian  lines  of  original  input  spacing 

IXZSYM  symmetry  option 

Ca 1  led  by 
NEW RAD 

A  - 2  9  Subroutine  NULYT 

Subroutine  NULYT  ..elects  a  subset  of  constant  u-velocity  and 
wing  and  pylon  thickness  panel  layout  data  to  be  saved  as  an  input 
iata  set  for  Program  II.  The  subroutine  also  calculates  and 
saves  the  net  strengths  associated  with  each  panel  corner.  A  list¬ 
ing  of  the  subroutine  is  presented  in  Figure  A-l(w),  and  a  flow  chart 
in  Figure  A  -  8 . 

The  datr  stored  in  the  new  arrays  describe  a  orid  of  corner 
points  representing  the  panel  layout,  thus  eliminating  some 
iuplica'ion.  For  example,  a  corner  which  is  the  junction  of  four 


panels  is  specified  four  times  in  the  panel  coordinate  arrays, 
only  once  in  the  corner  system.  Similarly,  the  slope  of  the 
boundary  of  two  panels  adjacent  in  a  chordwise  row  is  saved  only 
once  in  the  new  arrays.  However,  whenever  a  break  in  sweep  or 
dihedral  angle  occurs,  an  extra  chordwise  row  of  points  must  be 
saved  to  accommodate  the  velocity  calculations  to  be  performed 
in  Program  II. 

At  the  beginning  of  the  subroutine  indices  needed  to  describe 
the  new  arrays  are  calc  „ated.  Next,  as  the  flow  chart  indicates, 
data  is  stored  for  constant  u-velocity  panels  on  the  wing.  The 
sketch  below  illustrates  the  corner  numbering  convention  as  employed 
for  a  wing  with  one  leading-edge  sweep  break  and  six  panels. 

Panel  numbers  are  in  the  center  of  the  panels. 


For  each  chordwise  row  of  points,  the  variables  YPT, ZPT, SPHI,CPHI 
are  stored  only  once.  Variables  XPT  and  SWP  are  stored  for  all 
points.  The  slope  and  dihedral  values  used  correspond  to  panels 
to  the  left  of  the  points,  except  for  an  extra  row  of  points  which 
occurs  at  a  break  in  sweep  or  dihedral.  For  this  type  of  chord- 
wise  row  (points  7,8,9  in  the  above  sketch),  SPHI,  CHPI,  and 
SWP  correspond  to  panels  to  the  right. 

If  a  pylon  is  present,  the  subroutine  next  stores  the  corner 
data  for  pylon  panels.  These  data  are  stored  in  the  same  arrays 


following  the  wing  corner  data.  Similarly,  if  a  fuselage  is 
present,  the  corner  data  for  the  fuselage  panels  is  saved  in 
arrays  following  the  wing  and  pylon  data.  The  correspondence 
between  body  interference  panels  and  corner  point  numbering  is 
shown  in  th'  ketch  below  in  which  NBDCR1=3,  NBDCR2=2.  One 
should  nett  l.  =»t  an  extra  row  of  points  is  saved  for  each 
chordwise  row  of  panels. 
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The  subroutine  next  stores  the  data  for  wing  and  pylon  source 
panels,  using  the  same  procedure  as  for  wing  and  pylon  constant 
u-velocity  panels.  However,  the  source  panel  data  are  stored  in 
separate  arrays. 

The  last  section  of  the  subroutine  calculates,  by  means  of 
subroutine  NETCLC,  the  net  strength  associated  with  each  corner 
point.  Finally  at  the  end  of  the  routine,  a  factor  of  tt  which 
had  been  removed  from  the  source  panel  thickness  slopes  for 
boundary  condition  calculations,  is  restored. 

Subroutine  references 


NETCLC 


Called  by 
LDCALC 


A-30  Subroutine  PANEL 

Subroutine  PANEL'S  purpose  is  to  calculate  direction  cosines 
of  the  normal  vector  and  the  centroid,  area,  and  inclination 
angles  of  an  arbitrary  quadrilateral  panel.  It  is  called  to 
compute  the  geometric  properties  of  each  source  panel  on  the 
noncircular  body.  A  listing  of  the  routine  is  presented  in 
Figure  A-l(y)  of  this  report. 

The  four  corners  of  the  panel  are  numbered  in  a  clockwise 
direction.  A  diagonal  vector  connects  points  1  and  3,  and 
a  diagonal  vector  T^  connects  points  2  and  4.  The  normal  vector 
N  is  obtained  by  taking  the  cross  product  of  these  diagonal 
vectors,  and  the  direction  cosines  determined  by  calculating 
the  projections  of  this  vector  in  the  body  reference  coordinate 
system.  The  plane  of  the  panel  is  defined  to  be  perpendicular 
to  the  normal  vector  and  to  pass  through  a  point  whose  coordi¬ 
nates  are  the  averages  of  the  coordinates  of  the  four  input 
points.  The  input  points  are  then  projected  into  the  plane  of 
the  panel,  and  transformed  to  the  reference  coordinate  system. 

A  new  panel  coordinate  system  £ , n  is  introduced  with  the 
average  point  of  the  panel  as  origin.  The  coordinates  of  the 
centroid  and  the  panel  area  are  calculated  in  this  new  system, 
and  the  centroid  transformed  to  the  reference  system.  Two 
angles  are  used  to  define  the  panel  orientation.  The  incidence 
6  is  the  angle  between  the  x  axis  and  the  line  of  intersection 
with  the  panel  of  a  plane  passing  through  the  x  axis  and  per¬ 
pendicular  to  the  panel.  The  inclination  0  is  the  angle  between 
the  y  axis  and  the  line  of  intersection  of  the  panel  with  the 
yz  plane.  These  two  angles  are  calculated  in  terms  of  the 
direction  cosines  of  the  normal  vector. 


46 


The  descriptions  of  the  parameters  in  the  argument  list  follow 

J  panel  row  number 

K  panel  column  number 

NP  panel  number 

XYZ(NP,1),  array  containing  vectors  of  control  points,  XPT, 

XYZ(NP,6j  YPT,ZPT,  angles,  THET  and  DELTA,  and  area  for  NPth 
panel 

NBODY  total  number  of  panels 

XC , YC , ZC  array  containing  coordinates  of  corner  points  of 

panels 

KFUSOR  number  of  axial  stations  in  segment 

Subroutine  references 
INLTST,  INLXYZ 

Called  by 
B0DPN2 

A-31  Subroutine  PANVEL 

Subroutine  PANVEL  organizes  the  calls  to  SORPAN  for  the  calcu¬ 
lation  of  the  influence  of  a  panel  at  a  field  point.  This 
routine  computes  the  transformations  and  rotations  necessary  to 
calculate  the  influence  of  an  arbitrarily  oriented  panel  at 
either  a  field  point  or  at  the  control  point  of  another 
arbitrarily  oriented  panel.  A  listing  of  the  routine  is 
presented  in  Figure  A-l(z)  of  this  report. 

The  routine  first  computes  ...  ombined  transformations  to 
rotate  a  panel  or  its  image  into  the  local  panel  system  and  back. 
If  the  same  transformation  is  to  be  used  for  a  number  of  calcu¬ 
lations  the  logical  variable,  LZERO,  is  used  to  skip  around  this 
code  in  subsequent  calls.  The  field  point  in  question  is  rotated 
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into  the  panel  coordinate  system  and  a  call  made  to  SORPAN  to 
compute  the  influence  coefficient.  If  the  body  is  symmetric 
(IXZSYM=0)  the  above  calculations  are  repeated  for  the  symmetric 
panel  on  the  opposing  side  of  the  X-Z  plane.  The  influence  is 
summed  and  rotated  back  into  the  original  body  coordinate  system. 

The  descriptions  of  the  parameters  in  the  argument  list  follow: 

UB,VB,WB  component  velocity  influence  coefficients  in  body 

coordinate  system  at  field  point 
AN  resultant  normal  velocity  coefficient  of  panel  at 

a  field  point  or  at  another  arbitrarily  oriented 
panel 

IXZSYM  body  X-Z  plane  symmetry  option 

Subroutine  references 
SORPAN 

Called  by 

BODVEL,  FLDVL2 

A-32  Subroutine  PAS001 

Subroutine  PAS001  performs  the  L*U  decomposition  of  the 
positive  definite  matrix,  A.  This  routine  performs  no  pivot 
search  during  decomposition.  It  does  skip  unnecessary  calcula¬ 
tions  associated  with  off-diagonal  zeroes.  An  error  return  is 
provided  for  encountering  zero  values  on  the  diagonal.  The 
decomposition  procedure  is  equivalent  to 
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The  decomposed  matrix  is  stored  in  the  location  of  the  original 
matrix.  A  listing  of  the  routine  is  presented  in  Figure  A-l (aa)  of 
this  report. 

The  descriptions  of  the  parameters  in  the  argument  list  follow: 

positive  definite  matrix  with  nonzero  diagonal 
elements 

actual  size  of  matrix,  A 

diagnostic.  If  NER  greater  than  zero,  decomposition 
failed  because  | A (NER , NER) |  <  10-20 
dimensioned  size  of  A 

Called  by 
BODVEL 


A 

N 

NER 

ND 


A-33  Subroutine  PAS002 

Subroutine  PAS002  solves  the  system  of  equations  [L*U]  *  X  =  B  by 
forward  and  backward  substitution.  It  is  assumed  that  matrix.  A, 
has  been  decomposed  by  routine  PAS001  or  equivalent  so  that  A 
contains  [L*U],  A  listing  of  the  routine  is  presented  in  Figure 
A-l(aa)  of  this  report. 

The  descriptions  of  the  parameters  in  the  argument  list  follow: 

A  coefficient  matrix  containing  [L*U] 

B  matrix  containing  right-hand  sides  of  the  linear 

system  [L*U]  *  X  =  B.  Contains  X  on  output 
N  actual  size  of  [L*U]  matrix  contained  in  A 

NB  number  of  right-hand  side  vectors  contained  in 

the  first  NB  columns  of  B 
NDA  dimensioned  size  of  A 

NDB  dimensioned  size  of  B 


Called  by 
SMARCH 
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A-34  Subroutine  PLYOUT 


Subroutine  PLYOUT  reads  in  data  which  describe  the  geometric 
characteristics  of  the  pylon  and  calculates  quantities  which 
specify  the  pylon  constant  u-velocity  panels.  The  pylon  input 
variables  are  shown  in  Figure  8  of  Volume  II.  A  listing  of 
the  subroutine  is  presented  in  Figure  A-l (aa)  and  a  flow  chart  in 
Figure  A-9. 

The  first  part  of  the  subroutine  reads  in  and  prints  input 
items  43,  44,  and  45.  The  total  number,  MP,  of  pylon  constant 
u-velocity  panels  is  calculated.  Next,  the  subroutine  tests 
whether  the  pylon  leading-edge  and  trail ing-edge  sweep  angles 
are  constant  at  all  spanwise  z-stations.  If  breaks  in  sweep 
occur,  the  indicator  LVSPP  is  set  equal  to  one;  otherwise, 

LVSPP=0  and  the  quantity  SLPDIF,  the  difference  between  lead¬ 
ing-edge  and  trail ing-edge  slopes  for  a  chordwise  row,  is 
computed  outside  the  DO  loops.  The  zero  location  of  the  pylon 
root  chord  is  calculated,  and  the  variables  PLEX  and  CSIDE  are 
initialized  for  laying  out  the  first  chordwise  row  of  panels. 

The  remainder  of  the  subroutine  is  a  double  DO  loop  within 
which  panel  leading-edge  and  trail ing-edge  slopes,  corner  coordi¬ 
nates,  and  control  point  coordinates  are  calculated.  These 
quantities  for  the  pylon  are  stored  in  arrays  following  the  same 
quantities  calculated  by  subroutine  WLYOUT  for  the  NPANLS  constant 
u-velocity  panels  on  the  wing. 

The  following  sketch  shows  the  numbering  convention  associated 
with  the  pylon  panels  for  a  typical,  pylon-under-wing  layout  with 
four  chordwise  and  two  spanwise  panels. 


50 


In  the  subroutine  notation,  the  panel  corners  labeled  left  are 
those  farthest  from  the  pylon  root  chord. 

Called  by 
LDCALC 


A-3 5  Subroutine  PYTHIN 

Subroutine  PYTHIN  reads  in  the  pylon  thickness  data,  input 
items  46,  47  and  48.  A  listing  of  the  subroutine  is  presented 
in  Figure  A-l (bb)  of  this  report. 

The  number  of  panels  in  a  chordwise  row,  NCPS,  and  the  number 

of  chordwise  rows,  MSPS,  are  first  read  in  along  with  an  index, 

NUNIP,  which  indicates  whether  the  thickness  distribution  is 

similar  at  all  spanwise  stations.  The  subroutine  then  reads 

input  item  47.  Next,  if  NUNIP=1,  the  values  of  tan  6^  are  read 

in  for  the  first  chordwise  row  and  then  the  values  of  tan  9  for 

P 

the  other  rows  are  set  equal  to  those  of  the  first  row.  If  the 

distribution  is  not  similar,  NUNIP=0,  the  values  of  tan  9  for 

P 

all  rows  are  read  in. 


Called  by 
LDCALC 


A-36  Subroutine  RACKIO 


Subroutine  RACKIO  reads  and  prints  the  data  which  describe 
and  locate  the  rack  and  calculates  the  rack  source  and  doublet 
distributions.  The  rack  model  is  limited  to  description  as  a 
circular  body.  A  listing  of  the  routine  is  presented  in  Figure 
A-l(cc)  and  a  flow  chart  of  subroutine  RACKIO  is  presented  in 
Figure  A-10  of  this  report. 

The  subroutine  first  reads  and  prints  the  dimensional  lengths 
and  locations  of  the  rack,  and  the  polynomial  coefficients  which 
describe  the  shape.  The  rack  is  then  located  in  the  wing  coordi¬ 
nate  system.  Both  the  line  source  and  doublet  distributions  for 
the  rack  are  calculated  through  a  call  to  BDYGEN.  The  rack 
shoulder  is  located  from  the  first  zero  slope.  Lastly,  the  rack 
shock  wave  shape  is  computed  in  a  call  to  SHKSHP. 

Subroutine  references 
BDYGEN,  SHKSHP 

Called  by 
LDCALC 


A-37  Subroutine  SHAPE 

The  purpose  of  this  subroutine  is  to  calculate  the  body  radius 
and  surface  slope  at  a  specified  axial  station.  The  body  shape  is 
specified  by  a  series  of  polynomials  of  the  form  of  Equation  (1) 
of  Volume  II.  A  flow  chart  of  subroutine  SHAPE  is  presented  in 
Figure  A-ll  and  a  listing  of  the  subroutine  in  Figure  A-l (cc)  . 

The  calculation  performed  by  this  subroutine  consists  of  two 
steps.  The  first  step  is  to  determine  which  of  the  NS  polynomials 
describes  the  shape  at  the  value  of  X  where  the  radius  and  surface 


slope  are  required.  Once  this  is  determined,  the  appropriate  set 
of  coefficients  is  used  in  Equation  (1)  of  Volume  II  to  determine 
t/l.  The  value  of  dr/dx  is  found  by  differentiating  Equation  (1). 


It  should  be  noted  that  r/£  and  dr/dx  are  calculated  using  the 
coefficients  of  the  NSth  polynomial  if  x/l  is  greater  than  XE(NS). 

The  quantities  in  the  parameter  list  are: 

X  value  of  x/l  at  which  radius  and  surface  slope 

are  to  be  calculated 

NS  number  of  polynomials  describing  body  shape 

XE  array  containing  values  of  x/l  for  the  end  points 

of  the  NS  polynomials 

C  array  containing  the  coefficients  of  the  NS 

polynomials 

R  calculated  value  of  r/£  at  x/l  =  X 

DRDX  calculated  value  of  dr/dx  at  x/l  =  X 

Called  by 
BDYGEN 

A-38  Subroutine  SHKSHP 

Subroutine  SHKSHP  calculates  the  nonlinear  shape  of  the  shock 
wave  produced  by  an  axisymmetric  body  at  zero  angle  of  attack. 

The  calculations  performed  here  are  based  on  the  equations  pre¬ 
sented  in  Section  4.1.1  of  Reference  3.  A  listing  of  this  routine 
is  presented  in  Figure  A-l(dd)  of  this  report. 
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The  computation  sequence  used  here  is  to  incrementally  compute 
the  change  in  B  with  distance  from  the  centerline  and  numerically 
integrate  the  shape  radially.  The  procedure  starts  at  the  nose 
and  initializes  uhe  slope  of  the  shape  to  the  maximum  value 
BTNOSE  prescribed  through  the  input.  A  new  radial  distance  is 
prescribed,  and  an  axial  search  for  the  maximum  radial  velocity 
performed.  With  that  velocity  the  local  Pr andtl-Meyer  turning 
angle,  local  Mach  number,  and  new  shock  slope,  Bv.,  are  computed. 
The  new  axial  location  of  the  shock  shape  is  computed  from  the 
average  B  over  the  interval  step 
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The  integration  procedure  is  repeated  until  the  calculated  r> 

v/ 

approaches  the  free  stream  value  as  (1  -  3,^/6^)  <  0.01  or  a 
maximum  number  of  intervals  of  50  is  reached  before  terminating. 


The  descriptions  of  the  parameters  in  the  argument  list  follow 


NC 

SX(J) 
SS  (J) 
BETA 
RMAX 

BTNOSE 

NSHK 

XSHK 

RSHK 


maximum  number  of  axial  steps  used  to  scan  for 

maximum  radial  velocity 

axial  station  of  Jth  line  source 

line  source  strength  at  Jth  axial  station 

/  M^ -  1 ,  value  of  free  stream 

maximum  body  radius;  used  to  determine  initial 
step  size,  Ar 

limiting  value  of  B  at  nose;  input 
number  of  values  computed  in  shock  shape  table 
array  containing  axial  location  of  shock  relative 
to  nose 

array  containing  radial  location  of  shock  relative 
to  centerline 


Called  by 

FUSEIO,  RAC  K 1 0 ,  STORIO 
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A-39  Subroutine  SMARCH 


Subroutine  SMARCH  solves  for  the  source  panel  strengths  using 
a  ring  by  ring  marching  technique.  A  listing  of  the  subroutine 
is  presented  in  Figure  A-l (dd)  of  this  report. 

SMARCH  is  organized  to  solve  the  general  source  panel  solution 
for  a  body  in  supersonic  flow: 

[A] Yb  =  VB  -  [AI]fi  (A-14 ) 

where  [A]  is  the  aerodynamic  influence  coefficient  matrix 
partitioned  into  blocks  of  the  coefficients  of  one  ring  on 
another;  Yg  are  the  panel  strengths;  are  the  boundary  con¬ 
dition  normal  velocities  in  the  absence  of  an  image  body;  and 
[Aj]  is  the  influence  coefficient  matrix  of  the  image  body  on  the 
real  body  control  points.  No  image  store  effects  will  be  con¬ 
sidered  in  any  solutions  in  the  first  program.  The  solution 
proceeds  in  blocks  of  equations,  with  only  those  blocks  of 
equations  on  or  below  the  diagonal  computed  and  saved  on  TAPE9 . 

The  first  block  corresponding  to  the  influence  of  a  ring  on  itself 
is  read  from  TAPE9  in  L*U  decomposed  form.  The  solution  for  that 
ring  is  computed  with  PAS002.  Subsequent  blocks  in  column  form 
are  read  and  multiplied  by  the  strengths  of  that  ring  and 
subtracted  from  the  boundary  conditions  for  following  rings. 

The  descriptions  of  the  parameters  in  the  argument  list  follow 
GB  array  of  panel  strengths,  yB 

VB  array  of  panel  boundary  conditions  destroyed 

during  solution 
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IA  index  of  starting  location  in  blank  common  of 

temporary  A  matrix 

IMAGE  logical  indicator  of  presence  of  image  body 

Subroutine  references 
I OREAD ,  PASO  0  2 

Called  by 
SOLVE 


A-4 0  Subroutine  SOLVE 

Subroutine  SOLVE  computes  the  source  panel  boundary  condition 
for  a  noncircular  body  in  the  presence  of  a  free  stream  angle  of 
attack  and  roll  angle  and  calls  SMARCH  for  the  solution  for  the 
panel  strengths.  A  listing  of  the  routine  is  presented  in  Figure 
A-l(ee)  of  this  report. 

The  routine  initializes  the  angle  of  attack  parameters  and 
allocates  temporary  storage  in  blank  common  for  the  arrays 
required  during  the  calculation  of  panel  strengths.  The  boundary 
conditions  are  defined  as  follows 

V  =  cosa  sin6  -  cos6 ( s i na„cos$ cose  +  sinosimt  sine )  (A-15) 

The  solution  is  then  computed  in  SMARCH.  A  copy  of  blank  common 
up  through  the  panel  strengths  is  then  saved  on  TAPE7 .  SOLVE  may 
be  called  for  multiple  angles  of  attack  with  each  set  of  strengths 
saved  sequentially  after  the  other. 

The  description  of  the  parameter  in  the  argument  list  follows: 

IALP  index  of  the  number  of  the  IALPth  boundary  condition 

to  be  computed  and  saved.  IALP  <  1 
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S  u  b  t  o  u  t  l  r  1  e  r  i  •  f  <  ■  r  e  n  .  ■  e  s 
IObRIT.  SMARCH 


Called  by 

WDYEiDY,  STORIO 


A -4  1  Sub  tout  ine  b OK PAN 

Subroutine  SGRPAN  computes  the  three  components  of  velocity 
induced  at  a  specified  control  point  by  a  constant  source  dis¬ 
tribution  on  a  quadrilateral  panel  having  longitudinal  taper 
and  inclined  at  an  angle  delta  to  the  free  stream  direction. 

This  version  has  been  specialized  for  only  supersonic  flow. 

This  routine  is  based  on  the  methods  and  equations  presented 
in  Reference  4.  A  listing  of  the  subroutine  is  presented  in 
Figure  A-l(ff)  of  this  report. 

The  descriptions  of  the  parameters  in  the  argument  list  follow: 

U PM , VPM , WPM  three  orthogonal  components  of  velocity  in  local 

panel  coordinates  induced  by  panel  with  control 
point  XJ , ZJ  at  field  point  XI,YI,ZI 

Called  by 
PANVEL 


A-42  Subroutine  SOURCE 

Subroutine  SOURCE  calculates  coefficients  used  in  the  deter¬ 
mination  of  the  line  source  strengths.  They  occur  as  terms  in 
Equations  (1-14)  and  (1-17)  of  Reference  1.  The  relation  of  the 
coefficients  to  perturbation  velocities,  u  /V  and  v  /V  , 
induced  by  a  number  of  line  sources  distributed  along  the  body 
centerline  is  shown  in  the  first  two  of  Equation  (1-30)  in  which 
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the  coefficients  occur  as  the  multipliers  of  K  .  A  listing  of  the 
subroutine  is  presented  in  Figure  A-l(gg)  of  this  report.  The  sub¬ 
routine  is  called  by  subroutine  BDYGEN  and  Section  A-4  should  be 
referred  to  for  further  details  concerning  the  source  strength 
calculations. 

For  a  specified  control  point,  x  ,r  ,  and  singularity  origin, 

B  B 

',  the  subroutine  calculates  quantities  U  and  V  according  to  the 
following  equations: 


At  the  beginning  of  the  subroutine  a  test  is  performed  to  deter¬ 
mine  if  the  control  point  is  ahead  of  the  Mach  cone  from  the 
source  origin.  If  so,  U  and  V  are  set  tc  zero  and  control  is 
returned  to  the  calling  program. 

The  definition  of  the  variable  in  the  argument  list  is: 

TX  ;  location  on  body  axis  of  source  origin; 

positive  measured  from  tip  of  nose 

Called  by 
BDYGEN 

A-43  Subroutine  STORIO 

Subroutine  STORIO  reads  and  prints  the  data  which  describe  and 
locate  the  multiple  stores  and  organizes  the  calculation  of  the 
panel  or  line  singularity  strengths  for  both  circular  and  elliptic 
bodies.  The  routine  may  read  up  to  a  total  number  of  seven  stores 
of  either  circular  or  elliptic  shape.  Several  stores  may  share 
the  data  for  the  same  shape.  However,  when  elliptic  stores  are 


considered,  only  two  elliptic  shapes  are  permitted  with  up  to 
seven  stores  of  that  shape.  A  listing  of  the  subroutine  is 
presented  in  Figure  A-l (gg)  of  this  report.  A  flow  chart  of  the 
routine  is  presented  in  Figure  A-12. 

The  subroutine  consists  of  three  sequences  of  DO  loops.  The 
first  reads  and  prints  NSTRS  cards  containing  the  store  number, 
shape,  length,  maximum  radius,  location  relative  to  wing,  and 
orientation.  The  second  loop  reads  and  prints  NS  HPT  sets  of 
shape  data.  If  the  shape  is  circular,  items  57  through  59  of 
the  input  are  read  for  each  of  the  shapes.  If  the  shape  is 
elliptic,  items  60  through  75  are  read  through  a  call  to  GEOM. 

The  source  panel  geometry  is  comp  ited  and  a  temporary  copy  of 
all  shape  geometry  data  saved  on  TAPE11.  The  location  of  the 
store  nose  in  the  wing  cordinate  system  is  then  computed  for 
each  of  the  stores. 

The  last  major  loop  is  a  double  DO  loop  over  the  number  of 
shapes  and  over  the  number  of  stores.  The  outer  loop  is  over 
the  number  of  shapes  to  facilitate  reading  the  panel  geometry 
saved  on  TAPE11.  In  the  inner  loop  on  the  number  of  stores,  the 
store  shape  number  is  compared  with  the  shape  under  considera¬ 
tion  and  the  appropriate  properties  computed.  For  circular 
shapes,  BDYGEN  is  called  to  compute  the  line  source  and  doublet 
distributions,  the  shoulder  location  and  shock  shape  are  then 
calculated.  For  elliptic  shapes,  VELCMP  is  called  to  generate 
the  influence  coefficient  matrix.  Two  calls  to  SOLVE  may  be 
made.  For  the  first  store  of  the  shape  the  panel  strengths  are 
computed  at  zero  degrees  angle  of  attack  and  the  nonlinear  shock 
shape  computed  in  BSHOCK.  For  each  of  the  stores  of  that  shape 
SOLVE  is  called  again  to  compute  the  panei  strengths  at  the  angle 
of  attack  and  roll  of  the  fixed  store.  A  copy  of  the  store 
indices,  geometry  arrays,  and  strengths  are  then  saved  on  TAPE10 
by  FRSTRT. 


1 


Subroutine  references 

BDYGEN ,  BSHOCK,  FRSTRT ,  GEOM ,  SHKSHP,  SOLVE,  VEDCMP 

Called  by 
LDC  A  LC 


A-4 4  Subroutine  SWNGIN 

Subroutine  SWNGIN  reads  in  data  required  to  describe  the 
geometric  character i sties  of  the  wing  and  to  lay  out  the  wing 
constant  u-velocity  panels.  A  listing  of  the  subroutine  is 
presented  in  Figure  A-l(ii)  of  this  report. 

The  first  part  of  the  subroutine  reads  in  and  prints  input 
items  35,  36,  and  37,  which  consist  of  wing  geometry  data  and 
quantities  used  to  locate  the  trapezoidal-shaped  elemental  panels. 
If  any  nonzero  dihedral  angle  is  input,  the  logical  variable 
ZDIHED  is  set  equal  to  FALSE.  Next,  the  wing  twist  and  camber 
distribution,  if  any,  is  read,  input  items  38  and  39.  Two 
indices,  NTAC  and  NUNI,  are  first  input.  If  NTAC=0  there  is  no 
twist  and  camber.  The  index  NUNI  indicates  whether  the  twist  and 
camber  distribution  is  similar  at  all  spanwise  stations.  If  it  is 
similar,  NUN  I  =  1 ,  the  values  of  tan  are  read  in  for  the  first 
row  and  then  the  values  of  tan  a ^  for  the  other  rows  are  set  equal 
to  those  of  the  first  row.  If  the  distribution  is  not  similar, 

NUN  I  =0 ,  the  values  of  tan  for  all  rows  are  read  in. 

Called  by 

LDC  ALC 


A-4  5  Subroutine  THKLYT 


Subroutine  THKLYT  calculates  quantities  which  characterize 
the  wing  and  pylon  thickness  source  panels.  These  quantities  are 


stored  in  arrays  in  each  of  which  the  wing  panel  variables  precede 
the  pylon  variables.  All  panel  coordinates  are  expressed  in  the 
wing  coordinate  system  which  is  shown  in  Figure  7  in  Volume  II. 

A  listing  of  the  subroutine  is  presented  in  Figure  A-l(ii)  and  a 
flow  chart  in  Figure  A-13. 

The  first  part  of  the  subroutine  calculates  the  layout  of  the 
wing  panels.  After  the  variables  CSIDEP,  SLPDIF,  ZLER,  and 
WLEX  are  initialized  for  the  first  chordwise  row  of  panels,  the 
remaining  calculations  are  performed  within  a  double  DO  loop. 

The  outer  loop  index,  I,  specifies  the  chordwise  row  and  the 
inner  loop  index,  K,  the  panel  location  in  the  Ith  row.  Within 
the  inner  loop  the  panel  leading-edge  and  trail ing-edge  slopes, 
the  sine  and  cosine  of  the  dihedral  angle,  and  the  corner  coordi¬ 
nates  are  calculated  and  stored.  The  panels  are  numbered  con¬ 
secutively  in  chordwise  rows  beginning  with  panel  number  one  of 
row  one  adjacent  to  the  wing  root  chord  at  the  leading  edge.  The 
sequence  proceeds  in  increasing  numbers  to  the  trailing  edge,  then 
back  to  the  leading  edge  for  the  second  chordwise  row.  The  process 
continues  until  the  last  panel,  numbered  MS,  is  located  adjacent 
to  the  wing  tip  at  the  trailing  edge. 

If  a  pylon  is  present  (NPY=1),  the  remainder  of  the  subroutine 
calculates  the  panel  leading-edge  and  trail ing-edge  slopes  and  the 
corner  coordinates  for  the  pylon  source  panels.  The  procedure 
used  differs  very  little  from  the  wing  panel  calculations.  Values 
of  YPL,  ZPL,  and  LVSPP,  previously  calculated  in  subroutine  PLYOUT, 
are  used.  Also,  the  initial  value  of  PLEX  depends  upon  the  pylon 
location;  the  index  IP  is  tested  for  this  purpose.  The  numbering 
convention  associated  with  the  pylon  thickness  panels  is  the  same 
as  that  used  for  pylon  constant  u-velocity  panels.  It  is  described 
in  Section  A-34.  In  describing  corner  coordinates  for  both  wing 
and  pylon  panels,  the  corners  closer  to  the  respective  root-chords 


are  designated  right  corners;  those  farther  from  the  root-chords 
are  designated  left  corners. 

Called  by 
LDCALC 


A-46  Subroutine  THKOUT 

Subroutine  THKOUT  prints  the  slopes  of  the  wing  and  pylon 
thickness  distributions  which  were  read  in  as  items  42  and  48 
of  the  input  data,  see  Section  3.2.2.  After  the  thickness 
slopes  are  printed,  they  are  divided  by  it  and  saved  in  a 
combined  array,  DZDX,  for  subsequent  velocity  calculations. 

A  listing  is  presented  in  Figure  A-l(jj)  of  this  report. 


Called  by 
LDCALC 


A-4 7  Subroutine  VELBD1 

Subroutine  VELBD1  is  used  to  obtain  the  single-panel  u- 
velocity  panel  influence  coefficients  which  occur  indirectly  in 
Equations  (24)  ,  (25),  and  (26)  of  Reference  3.  Section  A-13 
of  this  report  should  be  consulted  for  further  details  concerning 
these  coefficients.  A  listing  of  the  subroutine  is  presented  in 
Figure  A-l(kk),  and  a  flow  chart  in  Figure  A-14  of  this  report. 
The  subroutine  uses  quantities  calculated  in  subroutine  BLYOUT 
which  is  described  in  Section  A-5.  The  wing  coordinate  system  is 
shown  in  Figure  7  in  Volume  II. 

At  the  beginning  of  the  subroutine  the  logical  variable  PYPNL 
is  set  equal  to  FALSE,  the  velocity  totals  UP,  VP,  and  WP  are 
initialized  to  zero,  and  the  panel  leading-edge  and  tr ai 1 ing-edge 
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slopes,  EM,  are  defined  as  zero.  In  the  remainder  of  the  sub¬ 
routine  the  influence  functions  for  the  four  corners  of  the  Ith 
panel  are  calculated  and  superposed.  The  superposition  scheme  is 
described  in  Section  3.3.4  of  Reference  3  and  that  reference 
should  be  consulted  for  further  details.  In  the  corner  numbering 
convention  for  fuselage  panels  corner  one  is  the  left  front, 
corner  two  the  right  front,  corner  three  the  left  rear,  and  corner 
four  the  right  rear  panel  corner.  The  right  corners  are  clockwise 
from  the  left  when  viewed  from  the  rear. 

Next  in  the  subroutine  the  corner  influence  function  totals 
TU,  TV,  and  TW  are  initialized  to  zero  and  a  test  is  performed  to 
determine  if  the  point  at  which  velocities  are  to  be  calculated 
lies  ahead  of  the  panel  leading  edge.  If  so,  all  calculations 
for  the  panel  are  skipped.  If  the  point  is  not  ahead  of  the 
leading  edge,  a  test  is  performed  next  to  determine  in  which 
fuselage  quadrant  the  panel  lies.  As  the  flow  chart  indicates, 
two  similar  but  distinct  transformation  and  superposition 
procedures  are  followed,  depending  upon  the  panel  location.  In 
each  procedure,  the  field  point  is  first  located  relative  to 
corner  one  in  a  corner  coordinate  system  which  is  shown  in  the 
sketch  below.  Each  panel  corner  has  an  associated  coordinate 
system.  The  xc  axis,  not  shown  in  the  sketch,  is  positive  to  the 
rear . 


If  corner  one  is  in  the  upper  left  quadrant  the  sign  of  the  yc 
coordinate  of  the  field  point  is  reversed  and  subroutine  VEL01  is 
called  to  calculate  the  influence  of  the  corner  on  the  image  of 

the  point  with  respect  to  the  x  ,z  plane.  The  sign  of  V,  which 

c  c 

is  the  returned  influence  function  in  the  y^  direction,  is  then 
reversed.  The  functions  U,V,W  are  resolved  back  into  the  wing 
system  and  superposed  in  the  same  manner  as  that  used  for  wing 
panels  with  positive  sweep  (see  Figure  4,  Reference  1).  Panels  in 
the  lower  left  quadrant  are  treated  in  the  same  manner  as  wing 
panels  with  negative  sweep.  No  further  coordinate  change  is 
necessary  and  subroutine  VEL01  returns  U,V,  and  W  which  are 
resolved  into  the  wing  system  and  superposed. 

Next,  in  each  procedure,  the  influence  of  the  mirror  image 
of  corner  one  with  respect  to  the  aircraft  vertical  plane  of 
symmetry  is  calculated.  This  is  accomplished  by  the  equivalent 
method  of  calculating  the  direct  influence  of  corner  one  on  the 
field  point  image.  Following  the  call  to  subroutine  VEL01  the 
sign  of  V  is  reversed  in  the  superposition. 

Corner  one  calculations  are  repeated  in  a  similar  manner  for 
corners  two,  three,  and  four,  but  only  corner  one  is  detailed  in 
the  flow  chart.  After  superposition  of  the  four  corner  influence 
functions  for  the  panel  is  completed,  the  influence  coefficients 
for  the  fuselage  panel  at  the  given  control  point  are  returned  by 
the  subroutine. 

Subroutine  references 
VE  L01 

Called  by 
DPCOEF 


A-48  Subroutine  VELCAL 


Subroutine  VELCAL  calculates  perturbation  velocities  at  a 
given  field  point  due  to  the  fuselage  and  store  source  and 
doublet  distributions.  A  listing  of  the  subroutine  is  presented 
in  Figure  A-l (mm)  of  this  report.  The  fuselage  coordinate  system 
is  shown  in  Figure  5  in  Volume  II.  The  store  coordinate  system 
is  shown  in  the  following  sketch. 


The  coordinates  of  the  field  point  are  given  as  formal 
parameters  in  the  appropriate  body  coordinate  system.  The 
subroutine  first  transforms  these  coordinates  into  the  VELCAL 
system  by  changing  the  sign  of  X  and  Z,  and  then  into  the  polar 
coordinates  XFIELD,  RFIELD,  and  THETA. 

The  major  part  of  the  program  consists  of  a  DO  loop  within 
which  the  axial,  radial,  and  tangential  velocities  due  to  the  N 
sources  and  doublets  are  calculated  and  summed.  A  test  is  made 
to  determine  whether  the  field  point  is  ahead  of  the  Mach  cone 
from  the  Ith  source  origin,  TX(I).  If  so,  all  remaining  perturba¬ 
tion  velocities  are  equal  to  zero  and  no  further  calculations  are 
performed  within  the  loop.  At  the  end  of  the  subroutine,  these 
velocities  in  the  VELCAL  system  are  resolved  back  into  the  body 
coordinate  system. 
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The  variables  in  the  subroutine  argument  list  are: 


T 

TC 

TX 


N 

XP 

Y 

ZP 
U  1 

VI 

W1 

BODYL 

SUMK 

SUMKD 


Called  by 
DPRHS 


array  containing  the  source  strengths 
array  containing  the  doublet  strengths 
array  containing  the  x  locations  of  the  origins 
of  the  singularities;  positive,  measured  aft  from 
tip  of  nose 

number  of  line  sources  and  doublets 
x  coordinate  of  field  point  in  body  system 

y  coordinate  of  field  point  in  body  system 

z  coordinate  of  field  point  in  body  system 

u/V  perturbation  velocity  at  field  point; 
body  system 

v/V^  perturbation  velocity  at  field  point; 
body  system 

w/V  perturbation  velocity  at  field  point; 
body  system 

body  x-station  at  which  base  singularities 
originate,  feet 

strength  of  source  originating  at  body  base 
strength  of  doublet  originating  at  body  base 


A-49  Subroutine  VELCMP 

Subroutine  VELCMP  organizes  and  calls  for  the  calculation  of 
the  aerodynamic  influence  coefficient  matrices  for  the  body  source 
panels.  This  routine  optionally  reads  the  Mach  number,  angle  of 
attac  and  angle  of  roll.  It  initializes  all  Mach  number  param¬ 
eters,  rewinds  any  external  files,  and  allocates  temporary  storage 
locations  in  blank  common  to  hold  blocks  of  influence  coefficients. 
VELCMP  performs  the  loop  over  the  number  of  body  segments  calling 
BODVEL  to  compute  the  influence  coefficients  for  the  panels  within 
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a  segment  on  all  other  control  points.  A  listing  of  this  sub¬ 
routine  is  presented  in  Figure  A-l (nn)  of  this  report. 

The  parameter  in  the  argument  list  is: 

IREAD  option  index  whether  to  read  the  Mach  number  and 

angle  of  attack 

IREAD=0 ,  read  XMAC H, ALPHAC , PH IR 

IREAD=1,  no  input  read;  parameters  must  be  in 

common 

Subroutine  references 

BODVEL,  IOREAD ,  IOWRIT 

Called  by 

WDYBDY,  STORIO 


A-50  Subroutine  VEL01 

Subroutine  VELOl  calculates  the  aerodynamic  influence 
functions  of  a  semi-infinite  triangle  associated  with  a  con¬ 
stant  u-velocity  panel,  as  described  in  Section  II-2.1  of 
Appendix  II  of  Reference  1.  The  influence  functions  relate 
the  panel  singularity  strength  to  the  perturbation  velocities 
induced  by  the  triangle  at  a  given  point.  They  occur  as  the 
coefficients  of  l/ir(u  /V^)  in  Equations  (II-4)  and  (11-12)  of 
Reference  1.  A  listing  of  the  subroutine  is  presented  in 
Figure  A-l (nn) ,  and  a  flow  chart  in  Figure  A-15  of  this  report. 

The  coordinate  system  used  by  the  subroutine  is  shown  in  Figure  3 
of  Reference  1. 

At  the  beginning  of  the  subroutine  the  logical  variable  FELT 
is  initialized  to  TRUE  and  a  test  is  performed  to  determine  if  the 
field  point  is  located  ahead  of  the  influencing  triangle  (X  <  0)  . 
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If  so,  the  influence  functions  U,V,W  are  set  to  zero,  FELT  is  set 
to  FALSE,  and  control  is  returned  to  the  calling  program. 


Next,  the  variable  PYPNL  is  tested  and,  if  the  triangle  is  on 
the  pylon,  a  transformation  is  performed  which  rotates  the  triangle 
into  the  VEL01  x,y  plane.  After  the  calculation  of  the  logical 
variable  INSIDE  and  other  frequently  used  quantities,  the  remainder 
of  the  subroutine  consists  of  four  major  sections  in  which  the 
influence  function  terms,  FI,  F2,  F4 ,  F5,  and  F7,  are  calculated. 
Each  section  corresponds  to  a  condition  of  the  slope,  EML,  of  the 
leading  edge  associated  with  the  sem i- in f i n i te  triangle.  The  sub¬ 
routine  requires  that  EML  ;>  0  and  this  is  accounted  for  in  the 
VEL01  calling  programs.  The  four  leading-edge  conditions  are 
described  fully,  with  accompanying  sketches  in  Section  II-2.1 
of  Reference  1.  All  equation  numbers  mentioned  in  the  following 
paragraphs  are  from  Section  II-2.1  of  Reference  1. 

The  first  section  of  the  subroutine  corresponds  to  a  subsonic 
leading  edge,  BTSQ  <  EMLSQ.  Equation  ( 1 1  —  5 )  is  used  if  the  point 
is  inside  the  Mach  cone  from  the  origin,  INSIDE=TRUE.  If  not, 

U,V,  and  W  are  set  to  zero.  In  this  section  as  in  the  remaining 
ones,  discontinuities  in  some  of  the  equations  may  occur  for 
certain  field  point  locations.  In  such  cases  the  affected 
influence  function  is  set  to  zero.  The  quantities  YYEDGE  and 
TLRNC ,  as  well  as  Y  and  Z,  are  used  to  test  the  singularity 
locations . 

If  BTSQ=EMLSQ,  the  leading  edge  is  a  sonic  leading  edge. 

The  equations  used  are  the  same  as  for  the  subsonic  case  except 
for  the  function  F2,  which  is  given  by  Equation  ( 1 1  —  7 ) .  If  the 
point  lies  outside  the  Mach  cone  from  the  origin,  the  influence 
functions  equal  zero. 
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The  third  section  of  the  subroutine  is  used  if  the  triangle 
leading  edge  is  supersonic,  BTSQ  >  EMLSQ.  Equations  (II-5)  and 
(II-9)  calculate  the  terms  of  the  influence  functions  if  INSIDE= 
TRUE.  If  not,  a  second  test  is  performed  and  Equation  (11-11)  is 
used  if  the  point  is  inside  the  Mach  cone  whose  origin  is  on  the 
leading  edge  at  the  field  point  y  location,  otherwise  the  func¬ 
tions  U , V , W  are  set  to  zero. 

The  fourth  section  of  the  subroutine  is  executed  if  the  leading 
edge  is  unswept,  EML=0.  For  this  special  case  the  perturbation 
velocity  equations  are  given  by  (11-12).  If  INS IDE=TRUE ,  the 
influence  function  terms  are  given  by  ( 1 1  —  5 )  and  (11-13).  Outside 
the  Mach  cone  from  the  origin  but  inside  the  cone  from  the  leading 
edge  Equation  (11-14)  is  used,  otherwise,  U,V,W  are  set  to  zero. 

The  last  part  of  the  subroutine  calculates  the  functions  U,V, 

W  from  the  component  terms,  in  the  case  of  a  leading  edge  with 
positive  sweep,  using  Equation  ( 1 1  —  4 )  .  If  the  triangle  is 
located  on  the  pylon,  V  and  W  are  rotated  back  into  the  pylon 
u-velocity  panel  coordinate  system. 

Called  by 

VELBD1,  VELPP1,  VELWP1 

A-51  Subroutine  VEL0T1 

Subroutine  VEL0T1  calculates  the  aerodynamic  inrlu  .  a  func¬ 
tions  of  a  semi-infinite  triangle  associated  with  wing  .r  pylon 
thickness  panel,  as  described  in  Section  II-2.2  of  Appendix  II  of 
Reference  1.  The  influence  functions  relate  the  panel  source 
strength  to  the  perturbation  velocities  induced  by  the  triangle 
at  a  given  point.  They  occur  as  the  coefficients  of  l/7i(tan  0) 
in  Equations  (11-15)  and  (11-16)  of  Reference  1.  A  listing  of 
the  subroutine  is  presented  in  Figure  A-l (pp)  of  this  report.  The 


subroutine  is  very  similar  in  logic  to  subroutine  VEL01  which  is 
described  in  detail  in  Section  A-50  and  represented  by  a  flow 
chart  in  Figure  A-15. 

In  subroutine  VEL0T1,  three  component  terms,  FI,  F2,  and  F5, 
need  to  be  calculated  in  order  to  determine  the  influence  func¬ 
tions  UTH,  VTH,  and  WTH.  Referring  to  Sections  II-2.1  and  II-2.2 
of  Reference  1,  function  FI,  F2,  and  F5  are  specified  in  Equation 
( 1 1  —  5 )  for  the  case  of  a  subsonic  leading  edge,  in  Equations 
( 1 1  —  5 )  and  ( 1 1  —  7 )  for  a  sonic  leading  edge,  and  in  Equations 
( 1 1 — 5 ) ,  (II-9),  and  ( 1 1  —  1 1 )  for  a  supersonic  leading  edge.  For 
the  special  case  of  an  unswept  leading  edge  ( EML=0 ) ,  the  general 
perturbation  velocity  equations  are  given  by  (11-16).  If  the 
given  point  lies  inside  the  Mach  cone  from  the  origin  of  the 
triangle,  function  FI  is  given  by  Equation  (11-13),  function  F5 
by  Equation  ( 1 1  —  5 )  ,  and  function  F2  by  Equation  (11-17).  If  the 
point  lies  outside  the  Mach  cone  from  the  origin  but  inside  the 
Mach  cone  from  the  triangle  leading  edge  at  the  field  point  y 
location,  functions  FI,  F2,  and  F5  are  given  by  Equation  (11-16). 

In  all  leading  edge  cases,  the  function  FI,  F2,  and  F5  are  singular 
for  certain  field  point  locations.  When  this  occurs,  the  affected 
influence  function  is  .et  to  zero.  All  influence  functions  are  set 
equal  to  zero  if  the  point  lies  in  the  plane  of  the  sem i - in f i n i te 
tr  iang le  ( ZTH=0 )  . 

A  table  equating  the  algebraic  and  program  notation  for 
variables  in  common  for  subroutine  VEL0T1  is  presented  in  Appendix 
B-2  of  this  report.  Almost  all  of  the  notation  used  in  subroutine 
VE LOT  1  is  defined  in  this  table.  One  should  note  that  the 
influence  functions,  U,V,W  and  the  point  coordinates,  YS,  ZS,  in 
VEL01  are  named  UTH,  VTH,  WTH,  YTH,  and  ZTH,  respectively,  in 
subroutine  VEL0T1. 

Called  by 

VE  LOT  1 ,  VE  LPT 1 


A-5_ 


Subroutine  VELPP1 


Su'^outine  VELPP1  calculates  the  single  panel  inf1 nonce  coef¬ 
ficients  which  form  the  coefficients  of  1/: (u+/v)  in  Equation 
(12)  of  Reference  1.  A  listing  of  the  subroutine  is  presented 
in  Figure  A-l  (qq)  and  a  flow  chart  in  Figure  A-16.  The  quan¬ 
tities  XX,  YY,  ZZ  in  the  subroutine  argument  list  are  the  x,y,z 
coordinates  in  the  wing  coordinate  system.  The  wing  coordinate 
system  is  shown  in  Figure  7  of  Volume  II. 

At  the  beginning  of  the  subroutine,  the  logical  variable 
PYPNL  is  set  equal  to  TRUE,  indicating  to  subroutine  VEL01  that 
calculations  are  to  be  performed  for  a  pylon  panel.  The 
quantities  YDIR  and  YIMG  are  calculated  and  the  velocity  totals, 

UP,  VP,  WP,  are  initialized  to  zero.  In  the  remainder  of  the 
subroutine  the  influence  functions  for  the  four  corners  of  the 
panel  are  calculated  and  superposed.  The  superposition  scheme  is 
described  in  Section  3.3.4  of  Reference  3  and  that  reference 
should  be  consulted  for  further  details.  The  corner  numbering 
convention  for  pylon  panels  associates  corners  one  and  two  with 
the  leading  edge  left  and  right  corners,  respectively;  corners 
three  and  four  with  the  trailing  edge  left  and  right  corners, 
respectively.  The  left  corners  are  those  farther  from  the  pylon 
root  chord. 

At  the  beginning  of  the  calculation  the  influence  function  totals 
are  initialized  to  zero  and  the  leading-edge  slope,  EMI,  and 
t r a i 1 i ng-edge  slope,  EM2,  are  defined.  Next,  a  test  of  the  sign 
of  EMI  is  performed  and  two  distinct  transformation  and  super¬ 
position  procedures  are  followeu  depending  on  the  results  of  this 
test,  as  the  flow  chart  indicates.  In  each  procedure  a  test  is 
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performed  first  to  determine  if  the  field  point  lies  abend  of  the 
most  forward  leading-edge  corner;  if  so,  all  calculations  for  the 
Ith  panel  are  skipped.  If  the  point  is  not  ahead  of  the  leading 
edge,  the  point  is  located  relative  to  corner  one  in  a  corner 
coordinate  system  which  is  illustrated  in  the  sketch  below. 

Each  panel  corner  has  an  associated  coordinate  system. 


For  a  panel  with  swept  back  leading  edge,  EMI  2  0  the  super¬ 
position  scheme  is  the  same  as  for  a  wing  panel  with  positive 
sweep  (see  Figure  4,  Reference  1).  The  field  point  z  coordinate  is 
reversed  and  subroutine  VEL01  is  called  to  calculate  the  influence 
of  corner  one  on  the  image  of  the  point  with  respect  to  the  x  ,y 
plane.  The  sign  of  W  is  then  reversed.  The  functions  U,V,W  are 
resolved  back  into  the  wing  system  and  superposed.  If  the  panel 
leading  edge  is  swept  forward,  superposition  is  the  same  as  for 
wing  panels  with  negative  sweep.  The  sign  of  EMI  is  reversed  and 
subroutine  VEL01  returns  U,V,W  which  are  resolved  into  the  wing 
system  and  superposed.  It  should  be  noted  that  the  subroutine 
code  combines  superposition  and  tr ansfo rmation  steps  and  that  the 
final  sign  changes  of  U  and  W  are  made  at  the  end. 

The  next  calculations  are  omitted  if  the  pylon  is  located  under 
the  fuselage  centerline.  Otherwise,  the  influence  of  the  mirror 
image  of  corner  one  with  respect  to  the  vertical  plane  of  symmetry 


is  calculated  .  This  is  accompl  i  si:-*d  by  t  i.<-  equi  v  .* ;  *  ;.*  <«■■.  ..  j 

calculating  the  direct  influence-  of  ~oi  n.-t  one  ■  n  1. 1„-  .  !.i  it: 

field  point  and  then  reversing  the  :  in  of 

Corner  one  calculations  are  rep.  ■•>•;*  .  i  n  a  ,• 

corners  two,  three,  and  four,  but  only  cot  nor  .*;>  i  :  i  n  ■  :  .. 

detail  in  the  flow  chart .  After  supet  pos n  ;•>;.«  . :  •  i..-  :  :: 

influence  functions  for  the  panel  is  comp;*-..  >  m-  . i  ; \h.  r. 

coefficients  for  the  Ith  pylon  panel  and  t  i  1  .  -  n  •  :  ,  i  p,;.-  i : .  t 

are  returned  by  the  subroutine. 

Subroutine  references 
VEL01 

Called  by 
DPCOEF 

A-5  3  Subroutine  VELPT... 

Subroutine  VELPT1  calculates  perturl.a-.iv  r*  y--  -.nci  t::os  a  l  a  f  i  > 

point  due  to  the  pylon  thickness  dist.rinr  orb  i j  to  mecho  : 

described  in  Section  3.3  of  Reference  ■  .  Th  •  .  ,  . .  am  logic  is 

very  similar  to  that  of  subroutine  VELPP  i  cui  c-.,l  ates  '--c  1  o- 

ities  induced  by  the  pylon  constant  u-veloc;  i.y  panels.  Subruut  in 
VELPP1  has  been  described  in  detail  in  Section  A-52  and  is  repre¬ 
sented  by  a  flow  chart  in  Figure  A-16  of  this  report.  Only  those 
details,  therefore,  in  which  subroutine  VELPT1  differs  from  sub¬ 
routine  VELPP1  are  included  in  this  description.  A  listing  of 
the  subroutine  is  presented  in  Figure  A-l(ss)  .  The  quantities 
in  the  subroutine  argument  list  are  the  x,y,z  coordinates  of  the 
field  point  in  the  wing  coordinate  system. 


The  coordinates  of  the  field  point  are  given  as  formal  param¬ 
eters  in  the  wing  coordinate  system.  The  point  is  located  rela¬ 
tive  to  each  of  the  four  panel  corners  using  the  same  transforma¬ 
tion  scheme  as  in  VELPP1.  However,  the  corner  coordinate  arrays 
which  define  the  pylon  thickness  panels  and  which  have  been 
previously  calculated  by  subroutine  THKLYT  (see  Section  A-45) 
are  used  in  the  transfo rmations.  Subroutine  VEL0T1  is  called 
to  calculate  the  corner  influence  functions,  U,V,W,  which  are 
then  superposed  in  the  same  manner  as  in  VELPP1. 

The  subroutine  calculates  perturbation  velocities  only.  The 
influence  coefficients  are  multiplied  by  DZDX(I)  to  obtain 
perturbation  velocities  induced  by  the  Ith  panel  at  the  given 
point.  These  velocities  are  calculated  and  summed  for  all 
pylon  thickness  panels. 

Subroutine  references 
VE LOTI 

Called  by 
DPRHS 


A-5 4  Subroutine  VELWP1 

Subroutine  VELWP1  calculates  the  influence  coefficients  at  a 
specified  control  point  due  to  a  wing  constant  u-velocity  panel. 
The  influence  coefficients  occur  in  Equation  (27)  of  Reference  3 
as  the  coefficients  oi  /;;)  (u,/V  ).  A  listing  of  the  subroutine 
is  presented  in  Figure  A-l(tt) ,  and  a  flow  chart  in  Figure  A-17. 
The  quantities  in  the  subroutine  argument  list  are  the  x,y,z 
coordinates  of  the  field  point  in  the  wing  coordinate  system. 

Tae  wing  coordinate  system  is  shown  in  Figure  7  of  Volume  II. 


At  the  beginning  of  the  subroutine,  the  logical  variable  PYPNL 
is  set  equal  to  FALSE  indicating  to  subroutine  VEL01  that  calcula¬ 
tions  are  for  a  wing  rather  than  a  pylon  panel.  The  influence 
function  totals,  TU ,  TV,  TW ,  are  initialized  to  zero.  The  leading 
edge  slope,  EMI,  and  the  trailing  edge  slope,  EM2,  as  well  as  the 
dihedral  angle  sine  and  cosine,  are  defined.  The  remainder  of  the 
subroutine  consists  of  the  calculation  and  superposition  of  the 
influence  functions  for  the  four  corners  of  the  influencing  panel. 
The  superposition  scheme  is  described  in  Section  3.3.4  of  Refer¬ 
ence  1  and  that  reference  should  be  consulted  for  further  details. 
The  corner  numbering  convention  for  wing  panels  is  shown  in 
Figure  4  of  Reference  1. 

Two  distinct  transformation  and  superposition  procedures  are 
followed  depending  on  the  sign  of  EMI,  for  corners  one  and  two, 
and  EM2 ,  for  corners  three  and  four.  In  each  procedure  the  point 
is  first  located  relative  to  corner  one  in  the  coordinate  system 
associated  with  subroutine  VEL01,  illustrated  in  Figures  3  and  4 
of  Reference  1.  Dihedral  effects,  if  nonzero,  are  included  in  the 
coordinate  transformations.  Subroutine  VEL01  is  called  to  calcu¬ 
late  the  corner  influence  functions,  U,V,W,  which  are  resolved 
back  into  the  wing  system  and  superposed.  Superposition  and 
transformation  steps  are  combined  in  the  code  and  the  final  sign 
changes  of  U  and  W  are  made  at  the  end  of  the  subroutint  .  Next, 
in  each  procedure,  the  influence  of  the  image  of  corner  one  with 
respect  to  the  aircraft  vertical  plane  of  symmetry  is  calculated. 
This  is  accomplished  if  EMI  <  0  by  the  equivalent  method  of 
calculating  the  direct  influence  of  corner  one  on  the  field 
point  image  and  then  reversing  the  sign  of  V. 

Corner  one  calculations  are  repeated  in  a  similar  manner  for 
corners  two,  three,  and  four,  but  only  corner  one  is  detailed  in 
the  flow  chart. 
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Subroutine  references 
VEL01 

Called  by 
DPCOEF 


A-55  Subroutine  VELWT1 

Subroutine  VELWT1  calculates  perturbation  velocities  at  a 
given  field  point  due  to  the  wing  thickness  distribution,  accord¬ 
ing  to  methods  described  in  Section  3.3  of  Reference  3.  The 
program  logic  is  very  similar  to  that  of  subroutine  VELWP1,  which 
calculates  influence  functions  for  the  wing  constant  u-velocity 
panels.  Subroutine  VELWP1  has  been  described  in  detail  in  Section 
A-54  and  is  represented  by  a  flow  chart  in  Figure  A-17  of  this 
report.  Only  those  details,  therefore,  in  which  subroutine  VELWT1 
differs  from  subroutine  VELWP1  are  included  in  this  description. 

A  listing  of  the  subroutine  is  presented  in  Figure  A-l (ww) . 

The  quantities  in  the  subroutine  argument  list  are  the  x,y,z 
coordinates  of  the  field  point  in  the  wing  coordinate  system. 

The  coordinates  of  the  field  point  are  given  as  formal  param¬ 
eters  in  the  wing  coordinate  system.  The  point  is  located  rela¬ 
tive  to  each  of  the  four  corners  of  the  trapezoidal  shaped 
thickness  panels  using  the  same  transformation  schemes  as  in 
VELWP1.  However,  the  corner  coordinate  arrays  which  define  the 
wing  thickness  panels  and  which  have  been  previously  calculated 
by  subroutine  THKLYT  (see  Section  A-45)  are  used  in  the  trans¬ 
formations.  Subroutine  VEL0T1  is  called  ^o  calculate  the  corner 
influence  functions,  U,V,W,  which  are  then  superposed  in  the 
same  manner  as  in  VELWP1. 

Finally,  the  subroutine  calculates  perturbation  velocities 
rather  than  single  panel  influence  coefficients.  Thus,  after 
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the  calculation  of  the  influence  coefficients  for  the  Ith  panel 
is  completed,  the  influence  coefficients  are  multiplied  by 
DZDX(I)  to  obtain  the  perturbation  velocities  induced  by  the 
Ith  panel  at  the  given  field  point.  These  velocities  are  calcu¬ 
lated  and  summed  for  all  wing  thickness  panels. 

Subroutine  references 
VE LOTI 


Called  by 
DPRHS 


A-56  Subroutine  WDYBDY 

Subroutine  WDYBDY  reads  and  prints  the  noncircular  fuselage 
data  and  organizes  the  calling  sequence  for  the  body  panel  layout, 
the  computation  of  the  source  strengths,  and  the  generation  of  the 
nonlinear  shock  shape.  A  listing  of  the  routine  is  presented  in 
Figure  A-l(zz)  of  this  report. 

WDYBDY  is  used  to  generate  the  required  solutions  for  the 
noncircular  fuselage.  It  first  reads  and  prints  the  basic 
length,  width,  nose  shock  angle,  and  interference  shell  length  and 
number  of  panels.  The  noncircular  fuselage  geometry  is  read  and 
panels  laid  out  in  a  call  to  GEOM.  The  panel  cross  section  data 
to  be  used  for  the  fuselage  interference  shell  is  saved  by  YZBIP. 
The  influence  coefficient  matrices  and  the  solution  for  the 
shock  strength  are  generated  in  calls  to  VELCMP  and  SOLVE.  The 
nonlinear  shock  shape  is  generated  in  a  call  to  BSHOCK  and  all 
data  to  be  passed  to  Program  II  saved  on  TAPE10  by  FRSTRT. 

Subroutine  references 


BSHOCK,  FRSTRT,  GEOM,  INLSHK,  SOLVE,  VELCMP,  YZBIP 


Called  by 
LDCALC 


A-57  Subroutine  WITHIN 

Subroutine  WITHIN  reads  in  the  wing  thickness  data,  input  items 
40,  41  and  42.  A  listing  of  the  subroutine  is  presented  in 
Figure  A-l(aaa)  of  this  report. 

The  number  of  panels  in  a  chordwise  row,  NCWS,  and  the  number 
of  chordwise  rows,  MSWS,  are  first  read  in  along  with  an  index, 
NUNIS,  which  indicates  whether  the  thickness  distribution  is 
similar  at  all  spanwise  stations.  The  subroutine  then  reads 
input  item  41.  Next,  if  NUNIS=1,  the  values  of  tan  6  are  read 
in  for  the  first  chordwise  row  of  panels  and  then  the  values  of 
tan  o  for  the  other  rows  are  set  equal  to  those  of  the  first  row. 

If  the  thickness  distribution  is  not  similar,  NUNIP=0,  the  values 
of  tan  9  for  all  rows  are  read  in. 

After  the  thickness  slopes  are  read  in,  an  input  error  check 
is  performed.  If  any  of  the  tan  0  values  is  zero  or  negative  at 
the  wing  leading  edge,  an  error  message  is  printed  (see  Section  3.5, 
Volume  II)  and  the  program  halts. 

Called  by 

LCALC 


A-58  Subroutine  WLYOUT 

Subroutine  WLYOUT  calculates  quantities  which  characterize 
the  constant  u-velocity  panels  on  the  wing.  An  example  of  a 
wing  trapezoidal  panel  is  shown  in  Figure  2  of  Reference  1. 
Similar  routines  PLYOUT  and  BLYOUT  calculate  these  quantities 
for  the  pylon  and  fuselage  panels,  respectively.  The  arrangement 
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of  variables  in  any  single  coordinate  array  is  wing  panels  first, 
pylon  panels  second,  and  fuselage  panels  last.  All  coordinates 
are  in  the  wing  coordinate  system  which  is  shown  in  Figure  7 
of  Volume  II.  A  listing  of  the  subroutine  is  presented  in 
Figure  A-l (aaa) ,  and  a  flow  chart  in  Figure  A-18  of  this  report. 

As  indicated  by  the  flow  chart,  the  first  part  of  the  sub¬ 
routine  tests  whether  the  wing  leading-edge  and  tr ai 1 ing-edge 
sweep  angles  are  constant  at  all  spanwise  Y-stations.  If 
breaks  in  sweep  occur,  the  indicator  LVSWP  is  set  equal  to  one; 
otherwise,  LVSWP=0  and  the  quantity  SLPDIF,  the  difference 
between  leading-edge  and  trail ing-edge  slopes  for  a  chordwise 
row,  is  computed  outside  the  main  DO  loop. 

After  initializing  the  quantities  CSIDEP,  ZLER,  and  WLEX, 
the  remainder  of  the  subroutine  consists  of  a  double  DO  loop. 

The  outer  loop  index,  I,  controls  the  chordwise  row;  the  inner 
loop  index,  K,  specifies  the  panel  location  in  the  Ith  row. 

Within  the  inner  loop,  the  leading-edge  and  trai 1 ing-edge 
slopes,  the  sine  and  cosine  of  the  dihedral  angle,  the  corner 
coordinates,  and  the  control  point  coordinates  are  calculated 
and  stored  for  each  panel.  The  convention  used  in  numbering 
the  NPANLS  wing  constant  u-velocity  panels  is  the  same  as  that 
used  for  the  wing  thickness  panels.  Section  A-45  should  be 
referred  to  for  details.  In  the  subroutine  notation,  the 
right-hand  side  of  the  panel  is  the  one  closer  to  the  root-chord. 

Called  by 

LDCALC 


A-59  Subroutine  WRFILE 

Subroutine  WRFILE  writes  the  data  set  on  TAPE12  which  contains 
all  the  information  required  to  continue  the  supersonic  store 
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1 


separation  computations  in  Program  II.  A  listing  of  the  routine 

is  presented  in  Figure  A-l(Lbb)  of  this  report. 

WRFILE  writes  all  ur  cata  required  for  Program  II 
unformatted  on  TAPE12.  The  data  includes  the  list  of  headings, 
control  indices,  various  constants,  u-velocity  data,  thickness 
panel  data,  fuselage  geometry  and  strength  data,  rack  data,  pylon 
data,  and  data  for  multiple  stores.  When  saving  information 
required  for  the  noncircular  fuselage  or  elliptic  stores,  the 
appropriate  arrays  are  transferred  from  TAPE10  to  TAPE12  using 
routine  FRSTRT. 

or- toutine  re i e r e n c e s 
F'KSTRT 

Called  by 
LDCALC 


A-6C  Subroutine  YZBIP 

Subroutine  YZBIP  is  useo  to  organize  the  scan  of  the  source 
panel  geometry  to  define  the  cross  section  data  for  the  non¬ 
circular  fuselage  interference  shell.  YZBIP  performs  the  loop 
over  the  number  of  body  segments  in  organizing  the  data  for  a 
scan  of  the  axial  geometry.  The  scan  of  panel  geometries  within 
a  segment  are  performed  in  a  call  to  YZBIP2.  A  listing  of  this 
subroutine  is  presented  in  Figure  A-l (ccc)  of  this  report. 

Subroutine  references 
YZ  SIP 


.’ailed  bv 
WDYBDY 
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A-61  Subroutine  YZBIP2 


Subroutine  YZBIP2  scans  the  source  panel  geometry  of  a  single 
body  segment  and  performs  the  interpolation  in  the  data  to  define 
the  y-z  values  of  the  section  to  be  used  by  BLY0T2  to  define  the 
body  interference  shell.  The  routine  scans  the  arrays  of  axial 
stations  to  find  the  two  sections  which  bound  the  input  station, 
XBIP  (item  15).  A  linear  interpolation  between  the  Y-Z  values  of 
the  two  sections  is  then  performed.  This  table  of  Y-Z  values  is 
used  by  BYL0T2  to  specify  the  Y-Z  values  of  the  fuselage  inter¬ 
ference  shell  as  close  to  the  source  panel  geometry  as  possible. 

A  listing  of  this  routine  is  presented  in  Figure  A-l (ccc)  of 
this  report. 

The  descriptions  of  the  parameters  in  the  argument  list  follow: 

XB  array  of  axial  stations  of  section  data  of 

segment 

YB , ZB  arrays  of  Y-Z  coordinates  of  panel  corners 

KRAD  number  of  corners  in  the  polar  direction  about 

the  body 

NFUSOR  number  of  axial  stations  in  segment 

Called  by 
YZBIP 
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Figure  A-l.-  Listing  of  Program  I 
(Pages  84  through  138) 
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(a) 

Figure  A-3.-  Flow  chart  of  subroutine  BDYGEN 


(a) 

Figure  A-4.-  Flow  chart  of  subroutine  BLYOUT. 


144 


Figure  A-5.-  Flow  chart  of  subroutine  DPCOEF. 


(b) 

Figure  A-6.- 


Continued 


rrn 


CALCULATE 
SOURCE  AMO 
DOUBLET 
DISTRIBUTION 


FIND  SHOULDER 
LOCATION 


CALCULATE  FUSELAGE 
SHOCK  WAVE  SHAi’E 


Figure  A-7.-  Flow  chart  of  subroutine  FUSEIO 


(b) 

Figure  A-9.-  Concluded 


Figure  A-10.-  Flow  chart  of  subroutine  RACKIO 


K  =  0 


I 


K  =  K  +  1 
XL  =  XE(K) 
J  =  K 


Figure  A-ll.-  Flow  chart 
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Figure  A-12.-  Flow  chart  of  subroutine  STORIO 


(b) 

Figure  A-12.-  Concluded. 
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(b) 

Figure  A- 16.-  Concluded 


(b) 

Fiqure  A- 17.-  Concluded. 
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APPENDIX  B 


COMMON  BLOCK  DESCRIPTIONS  -  PROGRAM  I 

B-l  Introduction 

The  purpose  of  this  Appendix  is  to  provide  more  detailed 
information  on  the  variables  passed  between  routines  through 
common  blocks  in  Program  I.  This  appendix  will  present  the  tables 
equating  program  notation  to  the  algebraic  notation  and  variable 
descriptions.  The  common  blocks  are  arranged  alphabetically.  It 
is  followed  by  a  section  detailing  the  special  usage  of  blank 
common.  When  a  description  identifies  an  item  as  an  input, 
consult  Volume  II  of  this  report  for  further  definitions.  For 
variables  used  both  in  the  description  of  the  noncircular  fuselage 
and  elliptic  stores,  only  the  input  item  for  the  fuselage  will  be 
referenced. 

A  cross-reference  chart  showing  the  routines  and  the  common 
statements  contained  in  each  is  presented  in  Figure  B-l.  Across 
the  top  of  the  chart  are  the  subroutine  names  including  the  main 
program  LDCALC .  Down  the  side  of  the  page  are  the  common  names. 

The  last  one  is  blank  common. 

B-?  Description  of  Variables  in  Labeled  Commons 

Var.  Engr.  Symbol  and  Description 

COMMON  /BASESG/  FBASE , FSUMK , FSUMKD, RBASE , RSUMK , RSUMKD , SBASE ( 7 ) , 
SSUMK ( 7 ) , SSUMKD ( 7 ) 

FBASE  fuselage  x-station  at  which  base  singularities 

originate,  feet 

FSUMK  strength  of  source  originating  at  fuselage  base 

FSUMKD  strength  of  doublet  originating  at  fuselage  base 

RBASE  rack  x-station  at  which  base  singularities 

originate,  feet 

RSUMK  strength  of  source  originating  at  rack  base 

RSUMKD  strength  of  doublet  originating  at  rack  base 


17.1 


SB  ASK 


store  x-station  at  which  base  singularities 
originate ,  tee t 

SSUMK  strength  of  source  originating  at  store  base 

SSUMKD  strength  of  doublet  originating  at  store  base 


COMMON  /BGEOM/  XFUS (51) , ZFUS ( 51 ) , FUSARD (51), FUSBY (51) ,FUSAZ(51) , 
X J  (51)  , PHIK (33) 

XFUS  (I)  Xq  coordinate  of  Ith  station  used  to  define  body 

external  geometric  shape,  feet;  input  item  18 

ZFUS(I)  zq  coordinate  of  Ith  station  containing  cambered 

offset,  feet;  input  item  19 

2 

FUSARD (I)  cross  sectional  area  at  Ith  station,  ft  ;  input 

item  22 

FISBY(I)  bw ,  elliptic  horizontal  semi-axis  { y-di rection ) , 

feet;  input  item  23 

FUSAZ(I)  a2,  elliptic  vertical  semi-axis  { z-direction ) , 

feet;  input  item  24 

XU (J)  XQ  coordinate  of  Jth  station  used  to  define  body 

panel  corners,  feet;  input  item  32 

RillK(J)  1^,  polar  angle  defining  Jth  meridian  of  panel 

edges,  degrees;  input  item  31 


COMMON  /BINLET/  N  INLET ,  NINVEL,  NT  INI, ,  RVIVO ,  NINBLK  ,  BTINLT  ,  YCP I  , 
XINLT , YINLT , ZINLT , XINLTE , YINLTE , ZINLTE , JINLT ( 25 ) 

NINLET  number  of  open  inlet  panels;  input  item  14 

NINVEL  number  of  additional  panels  to  be  used  in  velocity 

calculations  for  inlet  panels;  input  item  14 

MTiNL  total  number  of  inlet  panels  to  be  used  in  a  given 

calculation 

RVIVO  inlet  mass  flow  ratio;  ratio  of  open  inlet  panel 

frontal  area  to  total  inlet  panel  area 

NiNBLK  number  of  blocked  inlet,  panels;  input  item  14 

associated  with  inlet  panels;  input  item  29 


1  72 


BTINLT 


YCPI 


y- location  of  inboard  most  edge  of  inlet  panels 


XINLT , 
YINLT, 
ZINLT 


coor'di  nates  of  outboard  lead  i  nq  edge  of  inlet  panels 
used  to  locate  center  of  inlet  shock  propagation 


ZINLTE, 
YIN!  TE , 
ZINLTE 


coordinates  of  outboard  trailing  edge  of  inlet 
panels  used  to  locate  lower  lip  of  inlet  and 
turning  point  of  shock 


JINLT (I) 


fuselage  panel  number  associated  with  Jth  inlet 
panel 


COMMON  /BIMSHK/  NIS , XISHLD , EALPI , XCLOSD , MAXSH I , NINL ( 8 ) ,PHINL(3) , 
YINL ( 8 )  , X I N  L ( 8  0 )  ,RINL(80) 


NIS 

XISHLD 

EALPI 

XCLOSD 
MAXSH I 

NINL ( I ) 
PHINL(I) 

YINL (I) 

XINL 

RINL 


number  of  inlet  shock  traverses;  input  item  28 

x-station  at  which  . 1 s  associated  with  inlet  return 
to  free  stream 

angle  of  attack  correction  factor  for  inlet 
shock  (  E ALPHA) 

x-^  -cation  of  leading  edge  of  blocked  inlet  panels 

maximum  number  of  points  in  inlet  shock  tables  per 
traverse  ( =MAXSHK) 

number  of  points  computed  for  Ith  shock  traverse 

angle  measured  from  z-axis  below  inlet  of  Ith 
inlet  shock  traverse,  positive  counterclockwise 

y-station  of  Ith  inlet  shock  traverse 

table  of  x-values  of  inlet  shock 

table  of  radial  values  of  inlet  shock 


COMMON  /BIPYZ/  XBIP , HYZBIP , YCB ( 3 ‘ ' , ZCB ( 3 3 ) 

XBIP  x-station  at  which  the  y,z  geometry  used  in  the 

definition  of  the  body  interference  shell  are 
obtained  for  the  noncircular  fuselage,  feet 

NYZBIP  number  of  y,z  values  used  to  define  panel  corners 

around  the  circumference  of  the  interference  shell; 
NYZBIP=KRAD  of  local  body  segment 
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YCB(I)  y-coordinates  of  noncircular  interference  shell, 

feet 

ZCB(I)  z-coordinates  of  noncircular  interference  shell, 

feet 


COMMON  /BLKPAN/  COST , SINT , XBTJ , YBTJ , ZBTJ , XC1 , YC1 , ZC1 , XPTI , YPTI , 
ZPTI , COST I , SI NT I , COSD, SIND, LZERO 

COST, SINT  cos('jj)  and  sin(Oj),  cosine  and  sine  of  polar  angle 

of  Jth  influencing  source  panel 

XBTJ  XPTj  coordinate  of  Jth  influencing  control  point, 

feet 

YBTJ  Y~Tj  coordinate  of  Jth  influencing  control  point, 

feet 

ZBTJ  ZPTj  coordinate  of  Jth  influencing  control  point, 

feet 

XC1  X53  coordinate  of  Jth  panel  reference  corner,  feet 

YC1  yB  coordinate  of  Jth  panel  reference  corner,  feet 

ZC1  z  coordinate  of  Jth  panel  reference  corner,  feet 

XPTI  XFT-^  coordinate  of  Ith  influenced  field  point,  feet 

YPTI  YFT^  coordinate  of  Ith  influenced  field  point,  feet 

ZPTI  ZFT^  coordinate  of  Ith  influenced  field  point,  feet 

COSTI ,  SINTI  cos(S^)  and  sin(0j_),  cosine  and  sine  of  polar  angle 

at  Ith  influenced  panel  or  field  point 

COSD, SIND  cos(6j_)  and  sin(6j_),  cosine  and  sine  of  panel 

incidence  angle  at  Ith  influenced  panel  or  field 
point 

LZERO  PANVEL  angle  calculation  option:  F=yes,  T=no 


COMMON  /BODCOM/  AMACH , TAND , CX , XCOR ( 4 ) , YCOR ( 4 ) , ZCOR ( 4 ) , XI , YI , XJ , 
Z J , BETAO , BETAL , SUBSON , SUPERS 

AMACH  Mach  number  used  in  source  panel  influence 

calculation 
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TAND 


CX 

XCOR(K) 

YCOR(K) 

ZCOR(K) 

XI 

YI 

ZI 

XJ 

ZJ 

3ETA0 

DETAL 

SUBSON 

SUPERS 


tan6j,  tangent  of  incidence  angle  of  Jth 
influencing  panel 

panel  chord  length,  feet 

x  of  Kth  corner  in  local  panel  system,  feet 
y  of  Kth  corner  in  local  panel  system,  feet 
z  of  Kth  corner  in  local  panel  system,  feet 
x  of  Ith  field  point  in  local  panel  system,  feet 
y  of  Ith  field  point  in  local  panel  system,  feet 
z  of  Ith  field  point  in  local  panel  system,  feet 
x  of  Jth  panel  control  point  in  panel  system,  feet 
z  of  Jth  panel  control  point  in  panel  system,  feet 

/m2-1 ,  free  stream  Mach  number  constant 

00 

/M|- 1 ,  local  Mach  number  constant 

subsonic  logical  indicator;  SUBSON  =  AMACH.LT.l 

supersonic  logical  indicator;  SUPERS  =  AMACH.GT.l 


COMMON  /BOPTNS/  JO , J2 , J6 , NFUS , NRADX ( 5 ) ,NF0RX(5) , J2TEST, I PRES , ISOLV , 
INLET, IPLOT (4) ,IPRT(5) , IUVW , XSTART , XWLE , REFA , REFD , REFL , REFX , 
CCTEST , ITMAX , BODL , I Z 1 ( 12 ) 


JO 
J2 
J  6 

NFUS 
NRADX  ( I ) 

NFORX (I) 
J2TEST 

IPRES 


reference  area  indicator;  see  input  item  16 

body  type  indicator;  see  input  item  16 

body  camber  indicator;  see  input  item  16 

number  of  body  segments;  see  input  item  16 

number  of  points  used  to  define  section  of  Ith 
body  segment 

number  of  axial  station  on  Ith  body  segment 

parameter  to  specify  body  camber  and  cross-section 
definition 

not  used 


ISOLVE 


not  used 


INLET 


logical  inlet  indicator:  true  =  inlet  panels  present 

false  =  no  inlet  panel  present 


i  Plot 

IPRT 

IUVW 

XSTART 


XWLE 


not  used 

optional  print  control  parameter,  see  input  item  14 

component  velocity  calculation  option;  see  input 
item  14 

x-station  at  which  pressure  integration  is  started, 
feet 

x-station  at  which  pressure  integration  is  ended, 
feed 


REFA 


REFD 


REEL 

REFX , REFZ 

CCTEST 

ITMAX 

IZ1 


body  reference  area,  ft2;  see  input  item  17 

body  reference  length  used  for  moment  normalization, 
feet;  see  input  item  30 

body  length,  feet;  see  input  item  30 

x,z  coordinates  of  moment  reference  point,  feet; 
see  input  item  30 

solution  convergence  control  criteria  (=0.0001) 
solution  maximum  number  of  iterations  (=20) 
dummy  array,  not  used 


COMMON  /BPGEOM/  SNT2  (200)  ,CST2  (200)  , ZLC (200)  , ZRC (200)  , THTI  (200)  , 
SANGW 


SNT2 ( J) 
CST2 ( J) 
ZLC  (I) 

ZRC (I) 


sin(ANGj)  associated  with  Jth  body  panel 

cos(ANGj)  associated  with  Jth  body  panel 

coordinate  of  left  edge  of  Ith  constant 
u-velocity  panel 

Zyj  coordinate  of  right  edge  of  Ith  constant 
u-velocity  panel 


THTI  (J)  (s’  -  THTj)  ;THTj  is  value  of  THT  associated 

with  Jth  panel 

sin(ANGW)  =  ZBWO/FRMAX 
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SANGW 


COMMON  /BSHOCK/  NSHK(IC)  ,PHTS(10)  .TIICTNUO)  ,  MAXSIIK  ,  NSHOCK  ,  DBFTA  , 
i:  a  i  PiiA ,  ex  to ,  cnc;; ,  (  <■ ,  vs  nr  '  j  ou; ,  r.shk  ( i  go ; 

NSHK  ( I )  number  of  pointy  used  to  represent  f  th  modi  i  ied 

shock  shape 

PHIS  (I)  polar  angle  at.  which  1th  modified  shock  is  corputi' 

degrees;  input  item  J3 

THHTN  (I )  nose  limited  shoe1;  ingle  of  initial  shock  shape 

at  I  til  polar  angle,  degrees 

MAXSIIK  maximum  number  of  points  in  Itii  shock  shape;  input 

1 1  em  i 4 

NSHOCK  number  of  modified  shock  shape  computed;  input 

i tern  L4 

DBETA  not  used 

E  ALP  HA  angle  or  attack  correction  to  siiock  shape;  input, 

item  IS 

CNUO , CNU2  not  used 

XSHLDR  xB  location  of  body  nose  shoulder,  feet;  input 

item  15 

SHK  dummy  array,  not  used 

XSHK , RSHK  arrays  containing  xB  and  rB  locations  of  NSHOCK 

sets  of  NSHK(I)  points  representing  the  modified 
nose  shock  shape,  feet 


COMMON  /CAMBER/  ALPHAL(200) 

ALPHAL(J)  tan (  ) ,  slope  of  wing  camuerline  at  the  constant 

u-velocity  panel  control  points;  input  item  39 

COMMON  /CONFIG/  NFU , NPY , NSTRS , LVSWP , NRACK 

NFU  fuselage  indicator;  input  item  4 

NPY  pylon  indicator;  input  item  4 

NSTRS  number  of  stores  indicator;  input  item  4 


LVSWP 


NRACK 


breaks  in  wing  sweep  indicator;  LVSWP=0,  no; 
LVSWP= 1 ,  yes 

rack  indicator;  input  item  4 


COMMON  /CONSTS/  PI , PI2 , DTOR, RTOD, FOURPI 


DTOR 


RTOD 


ir/1 80  0 
180°/tt 


FOURPI 


COMMON  /DIMENS/  NX , NR, KX , KR, NXNR, KXKR, NXKR, MAXNX , MAXKR , NATOT , 

NBODY , KFUS , KRADX ( 5 ) ,KFORX(5) ,IXC(5) , IYC(5) , IZC(5) , IXZSYM, 
NADIM, NXDIM , NG , IXPT , IYPT , I ZPT ,  ITH , IDEL , NTAP7 , I AR , IAN, IUB, 

IGB  (7)  , I VB , IU , IV, IW , I VA , IWA , ICP , IPHI , I YB , NAG , NAP , NAV , NAS , 
NASHK , NAFLD , I AO, IDO , ISKO , IYIM, I ZIM , ISVN , ISKP , NRING , I ROW (50) 

NFUS 

NX  sum  of  axial  geometry  stations  (=  1  NFORX(I)) 

NR  not  used 

NFUS 

KX  sum  of  axial  panel  geometry  stations  (  =  1  KFORX(I)) 

KR  not  used 

NFUS 

NXNR  sum  (=  E  NFORX ( I ) *NRADX ( I ) ) 

NFUS 

KXKR  sum  (=  E  KFORX ( I ) *  KRADX ( I )  ) 

NFUS 

NXKR  sum  (  =  E  NFORX ( I ) *KRADX ( I ) ) 


KXKR 


NXKR 


MAXNX 


MAXKR 


NATOT 


maximum  of  (NFORX (I) , 1=1 , NFUS) 
maximum  of  (KFORX (I) , 1=1 , NFUS) 
last  location  accessed  in  blank  common 


NBODY 


KFUS 


KRADX (I! 


number  of  body  panels 

number  of  body  segments  paneled  (=NFUS) 

number  of  meridian  lines  used  to  define  panel  edges 
on  Ith  body  segment;  input  item  27 
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KFORX (I) 

IXC  (I) 

IYC  (I) 

IZC  (I) 

IXZSYM 

NADIM 

NXDIM 

NG 

IXPT 

IYPT 

IZPT 

ITH 
IDEL 
NTAP  7 
IAR 

IAN 

IUB 

IGB ( I ALP ) 

IVB 


number  of  axial  stations  used  to  define  leading  and 


trail ing 
item  27 

edges  of 

panels 

on 

Ith  body 

segment ; 

input 

location 
Ith  body 

in  blank 
segment 

common 

of 

start 

of 

XC 

array 

for 

location 
Ith  body 

in  blank 
segment 

common 

of 

start 

of 

YC 

array 

for 

location 
Ith  body 

in  blank 
segment 

common 

of 

start 

of 

ZC 

array 

for 

XZ-plane  symmetry  option;  input  item  14 
dimensioned  length  of  A  array  in  blank  common 
maximum  allowable  number  of  axial  stations  (=51) 


last  location  in  blank  common  containing  source 
panel  geometry  arrays 


location 

XPT 

in 

blank 

common 

of 

start 

of 

control 

points , 

location 

YPT 

in 

blank 

common 

of 

start 

of 

control 

points , 

location 

ZPT 

in 

blank 

common 

of 

start 

of 

control 

points , 

location 

in 

blank 

common 

of 

start 

of 

array  THET 

location 

in 

blank 

common 

of 

start 

of 

array  DELTA 

number  of  variables  last  written  on  TAPE 7 


location  in  blank  common  of  start  of  panel  areas, 

AREA 

location  in  blank  common  of  start  of  temporary  array 
AN  containing  influence  coefficients 

location  in  blank  common  of  start  of  temporary  array 
UB  containing  U,V,W  influence  coefficients 

location  in  blank  common  of  start  of  IALPth  array 
containing  the  source  strength  solution,  GB 

location  in  blank  common  of  start  of  temporary  array, 
VB,  containing  normal  velocity  boundary  conditions 
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TU 

location 

in 

blank 

common 

of 

start 

of 

U  velocities 

L  V 

location 

in 

blank 

common 

of 

start 

of 

V  ve 1 oci ties 

I  W 

location 

in 

blank 

common 

of 

start 

of 

W  velocities 

i VA, IWA 

locat ion 

in 

blank 

common 

of 

s  t  art 

o  f 

add  it i on a  1 

temporary 

velocity  array;. 

]  ■,’!> 

location 

in 

blank 

common 

of 

start 

of 

pressure 

coef  f icient 

,  CP, 

array 

I  Pill 

locati on 

in 

blank 

common 

of 

temporary 

PHI  array 

L  Y  B 

location  in  blank  common  of  start 
and  ZB  of  temporary  cross  section 

of  coordinates  YB 
geometry  in  NEWRAD 

:,A; , 

maximum 

locations 

in 

blank 

common 

required 

in 

GEOM 

NAP 

not  used 

NAV 

maximum 

locat ions 

in 

blank 

common 

required 

in 

VELCMP 

mas 

maximum 

locations 

in 

blank 

common 

required 

in 

SOLVE 

NASHK 

maximum 

locations 

in 

blank 

common 

required 

in 

BSHOCK 

NAFLD 

maximum 

locations 

in 

blank 

common 

required 

in 

FLDVEL 

I AO  offset  location  in  blank  common  of  all  above  arrays 

when  multiple  configurations  are  simultaneously  in 
core 


IDO  offset  location  in  blank  common  of  ID  array  containing 

/DIMENS/  information  for  additional  source  panel 
configurations  in  core 

ISKO  offset  location  in  blank  common  of  array  containing 

information  in  common  /BSHOCK/  for  additional  source 
panel  configurations  in  core 

IYIM, IZIM  location  in  blank  common  of  start  of  image  store 

Y  and  Z  control  point  coordinates 

ISVN, ISKP  location  in  blank  common  of  start  of  temporary 

arrays,  SVN  and  LSKP 

NRING  number  of  rings  of  panels  on  body 

IROW(I)  number  of  panels  around  Ith  ring  of  panels 
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COMMON 

/EXCEL/  UEI 

(20U) 

, via  (200)  , WEI  ( 

2  0 «)  .  , 

( :  I  i  d  0  u  i 

UK  1(1) 

uwi; 

, /v . 

ut  tiu_. 

<■  -  Lt.h 

-■out 

t't<  1  j  A  )  1  ! ;  * 

VE 1(1; 

vwL, 

.  /V, 

a '  the  , 

-  1 1  n 

c  o  n  t 

i  o  1  :  / 1 

WE  1(1) 

wWj 

J-  / 

../V., 

at  the 

=  Ith 

control  point 

CIR(l) 

arra 

y  con 

taininy 

right- 

hand 

s i tie s  of 

equations 

COMMON  /FLOW/  ALFACR,  GAMF ,  FMACII ,  RIIO  ,  V1NF  ,  BETA  ,  BETASg  ,  FMCllSg 


ALFACR 

lf ' 

radians 

GAMF 

not 

used 

i’MACll 

M  ; 

•.O 

input  item  3 

RHO 

not 

used 

VINF 

not 

used 

BETA 

3  = 

/m^-T 

uo 

BETASQ 

,2 

FMCHSQ 

M2 

oo 

COMMON  /FSGEOM/  FRMAX , SSPAN , FLTHC , BODYPL 
FRMAX  fuselage  maximum  radius;  input  item  5 

SSPAN  wing  semispan,  feet;  input  item  35 

FLTHC  fuselage  length,  feet;  input  item  5 

BODYPL  body  interference  shell  length,  feet;  input  item  1 


COMMON  /FSHOCK/  NFSHK , FXSHK ( 50 ) , FRSHK ( 50 ) , FDRDX ( 10 1 ) 


NFSHK 


number  of  x  and  r  values  in  fuselage  nose  shock 
table  generated  from  lino  singularities 


FXSHK  (I)  , 
FRSHK (I) 


x  and  r  coordinates  of  Ith  circular  fuselage  nose 
shock  location  at  zero  angle  of  attack 


FDRDX (J) 


DR/DX  at  Jth  circular  fuselage  control  point 
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COMMON  /FSOR/  FXL  (101)  ,  FSS  (100)  ,  FDS  (100)  ,  NFSOR 


FXL  array  containing  the  x  positions  of  the  fuselaye 

sources;  positive,  measured  from  tip  of  nose 

FSS  array  containing  the  strengths  of  the  circular 

fuselage  source  distribution 

FDS  array  containing  the  strengths  of  the  circular 

fuselage  doublet  distribution 

NFSOR  number  of  circular  fuselage  sources  and  doublets; 

input  item  9 


COMMON  /HEAD/  TITLE 1 ( 20 ) , TITLE 2 (20) 

TITLE1  array  containing  hollerith  description  of  non¬ 

circular  body  external  geometry 

TITLE2  array  containing  hollerith  description  of  non¬ 

circular  body  paneling  distribution 


COMMON  /ICVEL/  UP , VP , WP , I I , IF , DELTP ( 200 ) 

UP  u/V^  perturbation  velocity  at  a  control  point  due 

to  pylon  or  wing  thickness;  also  u-velocity  panel 
u  influence  function  at  a  control  point 

VP  v/Vx  perturbation  velocity  at  a  control  point  due 

to  pylon  or  wing  thickness;  also  u-velocity  panel 
v  influence  function  at  a  control  point 

WP  w/Vrjo  perturbation  velocity  at  a  control  point  due 

to  pylon  or  wing  thickness;  also  u-velocity  panel 
w  influence  function  at  a  control  point 

II, IF  initial  and  final  values,  respectively,  of  thickness 

panel  index 

DELTP  1/tt  (u+/Voo)  ;  u-velocity  panel  strength 

COMMON  /INDEX/  NCW , MSW , MSWP , NPANLS , NCWB , NBDCR1 , NBDCR2 , NBD , NBIP , 

MP , NCP , MSP , NIP , N2 , N2P , NPTOT 


NCW 

number 

of 

wing 

panels 

chordwise , 

input  item  36 

MSW 

number 

of 

wing 

panels 

spanwise , 

input  item  36 

MSWP 

MSW+1 

NPANLS 


number  of  u-velocity  panels  on  left  wing  panel: 
NCW*MSW 


NCWB 


NBDCR1 

NBDCR2 

NBD 

NBIP 

MP 

NCP 

MSP 

NIP 

N2 

N2P 

NPTOT 


number  of  DIP  rings,  input  item  9 

number  of  panels  in  ring  above  wing,  input  item  9 

number  of  panels  in  ring  below  wing,  input  item  9 

total  number  of  fuselage  u-velocity  panels  in  a 
ring  on  the  half  body 

number  of  u-velocity  panels  on  fuselage; 

NBIP  =  NBD*NCWB 

number  of  u-velocity  panels  on  pylon; 

MP  =  NCP  *MSP 

number  of  pylon  panels  chordwise,  input  item  44 
number  of  pylon  panels  spanwise,  input  item  44 
NPANLS  +  1 
NPANLS  +  MP 
N2  +  1 

NPANLS  +  MP  +  NBIP 


COMMON  /NUINDX/  KW , KP , NRW , NRP , NRB , NPTW , NPTP , NPTB , NRTOT , NPTT , 
NRWS , NRPS , NPTWS , NPTPS , NRTOT S , NPTTS , NRWP , NPTWP , 

NCW1 , NCP 1 , NCWB 1 , NCWS 1 , NCP  S 1 

KW  number  of  breaks  in  sweep  and/or  dihedral  on  wing 

KP  number  of  breaks  in  sweep  on  pylon 

NRW  NRW  =  MSW  +  KW  +  1 

NRP  NRP  =  MSP  +  KP  +  1 

NRB  NRB  =  2  *  NBD 

NPTW  NPTW  =  NCW  +  NRW  +  1 

NPTP  NPTP  =  NCP  +  NRW  +  1 

NPTB  NPTB  =  NCWB  +  NRB  +  1 

NRTOT  NRTOT  =  NRW  +  NRP  +  NRB 

NPTT  NPTT  =  NPTW  +  NPTP  +  NPTB 
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SINA 


sin (a) 


SINB  sin (3) 


COMMON  /PYGEOM/  Z ( 20 ) , XPLE , YPL, CRP , HP , PSIPLE ( 20 ) , PSIPTE ( 20 ) , IP , 
SLLE , PSLPOF , CENTER , ZPL , LVSPP 

Z (K)  Zp  locations  of  pylon  u-velocity  panel  side  edges; 

input  item  45 

XPLE  location  of  pylon  root  chord  leading  edge;  input 

item  43 


YPL 


yw  location  of  the  pylon;  YPL=Y(IP)  of  input 
item  37 


CRP 

HP 

PSIPLE (K) 

PHIPTE (K) 

IP 

SLLE 

PSLPDF 

CENTER 

ZPL 

LVSPP 


pylon  root  chord;  input  item  43 
pylon  height;  input  item  43 

leading-edge  sweeps  of  pylon  u-velocity  panels; 
input  item  45 

trailing-edge  sweeps  of  pylon  u-velocity  panels; 
input  item  45 

index  of  yw  location  of  the  pylon;  input  item  43 

slope  of  pylon  leading  edge  for  Ith  chordwise  row 
of  pylon  panels 

difference  in  leading  and  trailing  edge  slopes 
PSLPDF  =  SLLE-SLTE  for  Ith  chordwise  row  of 
pylon  panels 

logical  test  for  centerline  pylon  (CENTER=YPL. EQ. 0 ) 

z  location  of  the  pylon;  ZPL=ZLC ( (IP-2) *NCW+1) 
w 

breaks  in  pylon  sweep  indicator,  LVSPP=0,  no; 
LVSPP=1 ,  yes 


COMMON  /RKGEOM/  RRMAX , RLTHC , XWROC , YWRO , ZERO , NRPOLY , RXEND ( 7 ) , 
RCOEF (7, 7) , XBRO, YBRO , ZBRO , RIBCR, SRIBCR, CRIBCR 

RRMAX  maximum  rack  radius;  input  item  49 

RLTHC  length  of  rack;  input  item  49 
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XWROC 


YWRO 

ZWRO 

NRPOLY 

RXEND(J) 

RCOEF (I, J) 

XBRO 

YBRO 

ZBRO 

RIBCR 

SRIBCR 

CRIBCR 


xw  location  of  rack  in  wing  coordinates,  feet 

y  location  of  rack  in  wing  coordinates,  feet 

zw  location  of  rack  in  wing  coordinates,  feet 

number  of  rack  shape  polynomials;  input  item  50 
rack  polynomial  endpoints;  input  item  51 
rack  polynomial  coefficients;  input  item  52 
location  of  rack  in  body  coordinates,  feet 
yB  location  of  rack  in  body  coordinates,  feet 

zB  location  of  rack  in  body  coordinates,  feet 

RIODTOR;  input  item  49  for  RIC 
sin (RIBCR) 
cos (RIBCR) 


COMMON  /RSHOCK/  NRSHK, RXSHK (50) , RRSHK (50) , RDRDX (101) 

NRSHK  number  of  x  and  r  values  in  rack  nose  shock  table 

generated  from  line  singularities 

RXSHK (I),  x  and  r  coordinates  of  Ith  rack  nose  shock  location 

RRSHK (I)  at  zero  angle  of  attack,  feet 


RDRDX (J)  DR/DX  at  Jth  rack  control  point 


COMMON  /RSOR/  RXL (101) , RSS (101) , RDS (101) ,NRSOR 

RXL  array  containing  the  x  positions  of  the  rack 

sources;  positive,  measured  from  tip  of  rack  nose 

RSS  array  containing  the  strengths  of  the  rack  source 

distribution 

RDS  array  containing  the  strengths  of  the  rack 

doublet  distribution 


NRSOR 


number  of  rack  sources  and  doublets;  input  item  53 


COMMON  /SHPDAT/  NFPOLY , FXEND ( 7 ) , FCOEF (7,7) 

NFPOLY  number  of  fuselage  polynomials;  input  item  6 

FXEND  x/f  of  fuselage  polynomial  endpoints;  input 

item  7 

FCOEF  coefficients  of  fuselage  polynomials;  input 

ite.ii  8 


COMMON  /SSHOCK/  NSSHK(7) ,SXSHK(50,7) , SRSHK (50,7) , SDRDX (101,7) 

NSSHK(J)  number  of  x  and  r  values  in  Jth  store  nose  shock 

table  generated  from  line  singularities 

SXSHK ( I , J ) ,  x  and  r  coordinates  of  Ith  circular  store  shock 

SRSHK ( I ,  J )  location  at  zero  angle  of  attack  for  Jth  store 

SDRDX (K,J)  DR/DX  at  Kth  circular  store  control  point  of  Jth 
store 


COMMON  /SSOR/  SXL ( 101 , 7 ) , SSS ( 1 00 , 7 ) , SDS ( 100 , 7 ) , NSSOR ( 7 ) 

SXL ( I , J)  Ith  x-position  of  Jth  circular  store  sources; 

positive  measured  from  tip  of  nose 

SSS ( I , J )  Ith  source  strength  of  distribution  for  Jth 

circular  store 

SDS(I,J)  Ith  doublet  strength  of  distribution  for  Jth 

circular  store 

NSSOR (J)  number  of  Jth  circular  store  sources  and  doublets; 

equal  to  input  item  56 

COMMON  /STGEOM/  SRMAX ( 7 ) , SLTHC ( 7 ) , XWSOC ( 7 ) , YWSO ( 7 ) , ZWSO ( 7 ) , 
NSPOLY (7) , SXEND (7,7) , SCOEF ( 7 , 7 , 7 ) ,XBSO(7) ,YBSO(7^ , 

ZBSO ( 7 ) , SIBCR ( 7 ) ,NUMSTR(7) ,SSIBCR(7) ,CSIBCR(7) , 

SPHI RR ( 7 ) , NSHAPE ( 7 ) , NSHPT, MSHAPE (7) 

SRMAX (J)  maximum  radius  of  Jth  store;  input  item  54 

SLTHC ( J)  length  of  Jth  store;  input  item  54 

XWSOC (J)  xw  coordinate  of  tip  of  nose  of  Jth  store 

YWSO(J)  y  coordinate  of  tip  of  nose  of  Jth  store 


i 

\ 
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ZWSO(J) 

NSPOLY(J) 

SXEND (I , J) 

SCOEF (I , K, J) 

XBSO(J) 

YBSO(J) 

ZBSO(J) 

SIBCR(J) 

NUMSTR(J) 

SSIBCR(J) 

CSIBCR(J) 

SPHIRR(J) 

NSHAPE(J) 

NSHPT 

MSHAPE(K) 


z„  coordinate  of  tip  of  nose  of  Jth  store 
w 

number  of  polynomials  used  to  specify  the  shape 
of  the  Jth  store 

endpoints  of  polynomial  sections  specifying 
shape  of  Jth  store 

Jth  coefficient  of  Ith  polynomial  describing  the 
shape  of  the  Jth  store 

xB  coordinate  of  tip  of  nose  of  Jth  store 
yB  coordinate  of  centerline  of  Jth  store 
zB  coordinate  of  tip  of  nose  of  Jth  store 
SIC ( J) *DTOR;  input  item  54  for  SIC(J) 
store  number;  input  item  54 
sin (SIBCR ( J) ) 
cos (SIBCR (J) ) 

initial  roll  angle  of  Jth  store  coordinate  axes; 
positive  right  wing  down,  radians 

store  shape  number,  input  item  54 

number  of  different  store  shapes;  input  item  55 

shape  number  of  Kth  store  shape;  input  item  56 


COMMON  /SRCE/  XFIELD, RFIELD , BSQ , X2 , U, V, VT 


XFIELD 

RFIELD 

BSQ 

X2 

U 

V 

VT 


center  of  line  source/doublet  segment,  feet 
circular  body  radius  at  XFIELD,  feet 
not  used 
not  used 

axial  velocity  due  to  a  source 
radial  velocity  due  to  a  source 
not  used 
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COMMON  /THKDAT/  NCWS , NCPS , MS , MPS , MSWS , MSPS , NTHP , XRFT (400) ,XRBT(400) , 
XLFT (400) , XLBT ( 400 ) ,YRCT(400) ,YLCT(400) ,ZRCT(400) ,ZLCT(400) , 
THETAL ( 400 ) ,THETPL(200) ,SLLET(400) ,SLTET(400) ,DZDX(400) , 

YS ( 20 ) , PSWSLE ( 20 ) ,PSWSTE(20) ,PHIS(20) , ZS ( 20 ) ,PSPSLE(20) , 

P3PSTE (20) , SNPHS { 400 ) ,CSPHS (400) 


NCWS 

NCPS 

MS 

MPS 

MSWS 

MSPS 

NTHP 
XRFT ( I ) 

XRBT(I) 

XLFT ( I ) 

XLBT ( I ) 

YRCT(I) 

YLCT(I) 

ZRCT(I) 

ZLCT(I) 

THETAL ( I ) 

THETPL(I) 

SLLET(I) 


number  of  wing  thickness  panels  in  a  chordwise  row; 
input  item  40 

number  of  pylon  thickness  panels  in  a  chordwise  row; 
input  item  46 

number  of  thickness  panels  on  wing 

number  of  thickness  panels  on  pylon 

number  of  wing  thickness  panels  in  a  spanwise  row; 
input  item  40 

number  of  pylon  thickness  panels  in  a  spanwise  row; 
input  item  46 

number  of  thickness  panels  on  wing  and  pylon 

xw  coordinate  of  right  front  corner  of  Ith  thickness 
panel 

xw  coordinate  of  right  rear  corner  of  Ith  thickness 
panel 

x.  coordinate  of  left  front  corner  of  Ith  thickness 
panel 

xw  coordinate  of  left  rear  corner  of  Ith  thickness 
panel 

yw  coordinate  of  right  side  of  Ith  thickness  panel 

y  coordinate  of  left  side  of  Ith  thickness  panel 
w 

zw  coordinate  of  right  side  of  Ith  thickness  panel 

z.  coordinate  of  left  side  of  Ith  thickness  panel 
w 

wing  thickness  slopes;  input  item  42 
pylon  thickness  slopes;  input  item  48 
slope  of  leading  edge  of  Ith  thickness  panel 
slope  of  trailing  edge  of  Ith  thickness  panel 


SLTET(I) 


DZDX(I)  dz/dx  of  Ith  thickness  panel 

YS ( J)  yw  locations  of  wing  thickness  panel  side  edges; 

input  item  41 

PSWSLE(J)  leading-edge  sweeps  of  wing  thickness  panels; 

input  item  41 

:  SWSTE (J)  trailing-edge  sweeps  of  wing  thickness  panels; 

input  item  41 

PHIS ( J)  wing  thickness  panel  dihedrals;  input  item  41 

ZS(J)  Zp  locations  of  pylon  thickness  panel  edges; 

input  item  47 

PSPSLE(J)  leading-edge  sweeps  of  pylon  thickness  panels; 

input  item  47 

PSPSTE(J)  trailing-edge  sweeps  of  pylon  thickness  panels; 

input  item  47 

SNPHS(J)  sin (PHIS (J)) 

CSPHS(J)  cos (PHIS (J)) 

COMMON  /THVARG/  X,  YV,  ZV,  D,  W,WV,EML,PRT* 

X , YV, ZV  xfyfz  coordinates  of  field  point  in  thickness 

panel  coordinate  system 

U,W,WV  u,v,w  velocity  components  in  thickness  panel 

coordinate  system 

EML  slope  of  thickness  panel  leading  or  trailing  edge 

PRT  not  used 

COMMON  /VELARG/  X , YV, ZV, U, W, WV, EM, TLRNC , TIPY , PYPNL* 

X, YV, ZV  x,y,z  coordinates  of  field  point  in  u-velocity 

panel  coordinate  system 

U,W,WV  u,v,w  velocity  components  in  u-velocity  panel 

coordinate  system 


★ 

Variable  names  may  differ  from  routine  to  routine  but  the 
definitions  are  unchanged. 
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EM  slope  of  u-velocity  panel  leading  or  trailing  edge 

TLRNC  error  tolerance 

TIPY  exposed  wing  span 

PYPNL  logical  pylon  indicator 


COMMON  /WGEOM/  XBWCC, ZBWO, CRW, SLPWLE, SLPWTE , PSIWLE (20) , 

PSIWTE ( 20 ) ,  Y (20) ,PHID(20) ,ZDIHED,WICR 

XBWOC  xB  coordinate  of  wing  root  chord  leading  edge; 

input  item  34 

ZBWO  zB  coordinate  of  wing  root  chord  leading  edge; 

input  item  34 

CRW  wing  root  chord;  input  item  35 

SLPWLE  tan (PSIWLE (I) ) 

SLPWTE  tan (PSIWTE (I) ) 

PSIWLE (I)  leading-edge  sweeps  of  wing  u-velocity  panels; 

input  item  37 

PSIWTE (I)  trailing-edge  sweeps  of  wing  u-velocity  panels; 

input  item  37 

Y ( I )  yw  locations  of  wing  u-velocity  panel  side  edges; 

input  item  37 

PHID(I)  dihedral  angles  of  wing  u-velocity  panel  sections; 

input  item  37 

ZDIHED  logical  variable  indicating  whether  or  not  there 

is  wing  dihedral;  ZDIHED=TRUE,  no  dihedral; 

=FALSE ,  there  is  dihedral 

WICR  wing  incidence  angle,  radians;  computed  from  WIC, 

input  item  34 


COMMON  /WPGEOM/  XRF ( 200 ) , XRB ( 200 ) , XLF ( 200 ) , XLB ( 200 ) , YRC (200 ) , 

YLC (200) ,SCPT (200) ,YCPT(200) ,ZCPT(200) ,SWPPLE(200) , 

SWPPTE ( 200 )  ,SNPHI  (200)  ,CSPHI  (200) 

XRF (I)  xw  coordinate  of  right  rear  corner  of  Ith  constant 

u-velocity  panel 

XRB (I)  xw  coordinate  of  left  front  corner  of  Ith  constant 

u-velocity  panel 
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XLF(I) 


xw  coordinate  of  left  rear  corner  of  Ith  constant 
u-velocity  panel 


XLB(I)  xw  coordinate  of  right  front  corner  of  Ith  constant 

u-velocity  panel 

YRC(I)  yw  coordinate  of  left  edge  of  Ith  constant  u-velocity 

panel 

YLC(I)  yw  coordinate  of  right  edge  of  Ith  constant 

u-velocity  panel 

XCPT(I)  xw  coordinate  of  Ith  control  point 

YCPT(I)  yw  coordinate  of  Ith  control  point 

ZCPT(I)  zw  coordinate  of  Ith  control  point 

SWPPLE ( I )  sweep  of  leading  edge  of  Ith  constant  u-velocity 

panel  on  wing,  fuselage  or  pylon 

SWPPTE ( I )  sweep  of  trailing  edge  of  Ith  constant  u-velocity 

panel  on  wing,  fuselage  or  pylon 

SNPHI(I)  sine  of  dihedral  angle  of  Ith  constant  u-velocity 

panel  on  wing 

CSPHI(I)  cosine  of  dihedral  angle  of  Ith  constant  u-velocity 

panel  on  wing 

COMMON  /XSHOLD/  FXSHLD , RXSHLD , SXSHLD ( 7 ) 

FXSHLD  xD  location  of  circular  fuselage  shoulder,  feet 

D 

RXSHLD  xR  location  of  rack  shoulder,  feet 

SXSHLD (J)  xs  location  of  Jth  circular  store  shoulder,  feet 


B-3  Blank  Common 


The  requirement  for  handling  the  many  large  arrays  associated 
with  the  solution  for  u-velocity  and  source  panel  strengths  has 
necessitated  both  the  use  of  out  of  core  data  handling  and  storage 
and  the  setting  aside  of  a  scratch  storage  area  in  core  to  be  used 
for  more  than  one  purpose.  To  handle  the  latter  data  requirement 
blank  common  has  been  reserved  for  all  calculations  involving 
large  arrays.  The  program  flow  of  calculations  is  thus  arranged 
to  allow  variables  to  be  read  from  or  written  to  external  files 
as  needed.  The  following  describes  the  Program  I  sequence  of 
references  to  external  files  and  which  arrays  reside  in  blank 
common  at  each  point  in  the  program  flow.  The  information  residing 
in  each  of  the  external  files  is  described  later. 

In  Program  I  blank  common  is  used  for  three  purposes:  (1) 
storing  the  dynamically  dimensioned  arrays  associated  with  both 
the  noncircular  fuselage  and  the  elliptic  store;  (2)  storing  the 
wing-fuselage  u-velocity  aerodynamic  influence  coefficients  during 
the  solution  for  panel  strengths;  and  (3)  temporary  storage  of 
arrays  associated  with  the  summation  of  u-velocity  singularity 
strengths  at  panel  corners.  The  descriptions  which  follow  focus 
on  the  definition  of  quantities  during  these  phases.  The  non¬ 
circular  fuselage  and  elliptic  store  arrays  are  lumped  together 
because  they  share  common  code.  Only  the  dimensions  of  arrays 
will  vary  between  components.  The  arrangement  of  variables  in 
blank  common  during  these  calculations  may  change  dynamically  as 
the  solution  progresses  as  temporary  arrays  are  required  or  dis¬ 
carded.  The  flow  chart  in  Figure  A-2  of  Appendix  A  also  identifies 
several  points  in  the  program,  item  numbers,  where  external  files 
are  referenced.  The  comments  which  follow  are  keyed  to  the  usage 
of  blank  common  at  those  points  wherever  possible.  The  descrip¬ 
tions  of  the  use  of  blank  common  during  these  three  phases  follow. 

The  first  use  of  blank  common  is  to  store  the  arrays  of  panel 
properties  and  solutions  of  either  the  noncircular  fuselage  or  the 


elliptic  store.  This  spans  items  1  through  4  for  the  fuselage  and 
repeats  as  item  5  for  the  store  in  Figure  A-2.  The  first  routine 
using  blank  common  is  GF.OM  as  called  from  WDYBDY.  The  noncircular 
and  elliptic  store  both  share  the  same  analysis  routines  and  blank 
common  data  arrays.  In  Program  I,  the  results  of  the  first  are 
computed  and  all  information  required  for  Program  II  are  saved. 

The  solution  for  a  single  body  is  broken  into  four  parts.  These 
follow  the  progress  of  calculations  and  correspond  to  generation 
of  panel  geometry,  computation  of  panel  aerodynamic  influence 
coefficients,  solution  of  the  equations,  and  saving  of  appropriate 
arrays  required  for  Program  II.  Because  of  the  large  number  of 
possible  combinations  of  paneling  schemes,  all  array  storage 
allocation  m  blank  common  is  performed  dynamically.  That  is,  the 
array  length  is  computed  and  core  locations  set  aside  prior  to 
calculation  of  the  array  itself.  Arrays  are  typically  stacked  one 
behind  another  with  no  unused  locations  in  between.  No  explicit 
array  names  or  dimensions  are  defined  in  blank  common.  Only  the 
address  in  blank  common  of  the  first  element  in  the  array  is  saved 
in  labeled  common,  DIMENS. 

The  first  use  of  blank  common  within  the  calculation  of  the 
source  panel  geometry  (item  1)  is  for  temporary  storage  of  the 
arbitrary  Y,Z  input  of  the  external  fuselage  shape  as  read 
directly  from  input  cards.  They  are  followed  by  the  intermediate 
YB  and  ZB  values  computed  at  axial  stations,  XFUS ,  at  the  new 
meridional  angles,  PHIK,  of  the  revised  panel  spacing  layout. 

Blank  common  is  equivalently  dimensioned  to: 

COMMON  SFUS (NRAD , 2 , NFUSOR) , YB (NFUSOR, KRAD) , ZB (NFUSOR, KRAD) , . 

A(NADIM) 

where 

SFUS (NN , 1 , N )  y-station  in  local  body  coordinates  of  the 
NNth  meridian  at  the  Nth  axial  station, 

XFUS 

SFUS (NN , 2 , N)  z-station  in  local  body  coordinates  of  the 
NNth  meridian  at  the  Nth  axial  station, 

XFUS 


194 


YB(N,K) 


y-station  in  local  body  coordinates  of  the  Nth 
axial  station  for  the  Kth  revised  meridional 
angle 

ZB (N, K)  z-station  in  local  body  coordinates  of  the  Nth 

axial  station  for  the  Kth  revised  meridional 
angle 

If  more  than  one  segment  exists  YB  and  ZB  are  saved  on  TAPE8 
for  each  of  NFUS  segments. 

The  second  change  in  the  use  of  blank  common  as  identified  by 
item  2  in  Figure  A-3  is  to  store  the  arrays  containing  the  geometric 
properties  of  the  source  panels.  The  configuration  of  blank  common 
at  this  point  in  the  calculations  remains  the  same  while  con¬ 
tinuing  to  grow.  Arrays  which  are  to  be  saved  are  located  at  the 
beginning  with  scratch  space  at  the  rear.  From  the  end  of  the 
calculations  in  GEOM,  the  arrays  in  blank  common  contain  all  the 
panel  geometric  properties  required  for  calculation  of  influence 
coefficients  and  the  resolution  of  forces  and  moments.  Additional 
temporary  arrays  are  allocated  to  hold  the  ring-by-ring  coefficient 
arrays  in  VELCMP.  The  equivalently  dimensioned  arrays  in  blank 
common  look  like: 


COMMON  XPT(NBODY) , YPT (NBODY) ,ZPT(NBODY) ,  THET (NBODY) , 
DELTA ( NBODY ) , AREA (NBODY) ,XC(KX) ,YC(KXKR) , 
ZC(KXKR) , AN (MAXKR, MAXKR) , UB ( 3 , MAXKR, MAXKR) , . . . , 
A(NADIM) 

where 


XPT(I) 
YPT (I) 
ZPT(I) 
THET (I) 
DELTA ( I ) 
AREA (I) 
XC (L,M) 


x-station  of  control  point  of  Ith  panel 

y-station  of  control  point  of  Ith  panel 

z-station  of  control  point  of  Ith  panel 

inclination  angle  at  Ith  panel  control  point 
incidence  angle  of  Ith  panel 
surface  area  of  Ith  panel 

panel  corner  points  at  Lth  axial  station  of 
Mth  segment 
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YC(L,N,M) 


y  corner  point  at  Lth  axial  station  of  Nth 
meridional  angle  of  Mth  segment 


ZC (L , N , M)  z  corner  point  at  Lth  axial  station  of  Nth 

meridional  angle  of  Mth  segment 

AN  ( I , J)  temporary  array  for  aerodynamic  influence 

coefficients 


UB  ( 1 , 1 , J )  temporary  u-component  influence  coefficient 

array 

UB  ( 2 , 1 , J )  temporary  v-component  influence  coefficient 

array 


UB  ( 3 , 1 ,  J ) 


temporary  w-component  influence  coefficient 
array 


The  dimensions  and  indices  defining  the  starting  locations  of 
each  of  these  equivalent  arrays  are  found  in  common  DIMENS.  The 
starting  indices  may  typically  be  found  by  preceding  the  first  two 
or  three  letters  of  the  array  name  with  I. 

Additional  arrays  associated  with  the  solution  for  the  panel 
source  strengths  in  SOLVE  as  identified  by  item  3  in  Figure  A-2. 
Only  the  arrays  containing  the  panel  strengths  are  saved.  The 
strengths  for  up  to  seven  flow  conditions  may  be  retained  to  allow 
for  using  the  same  shape  data  for  more  than  one  store.  The 
equivalently  dimensioned  arrays  in  blank  common  looks  like: 

COMMON  XPT(NBODY) , YPT (NBODY ) , ZPT (NBODY ) , THET (NBODY) , 

DELTA (NBODY ) , AREA (NBODY) ,XC(KX) ,YC(KXKR) ,ZC(KXKR) , 

GB (NBODY, IALP) ,VB (NBODY) , AN (MAXKR, MAXKR) , . . . , 

A (NADIM) 

where  the  arrays  XPT  through  ZC  contain  the  same  information  as 
before  and 

GB(I,J)  strength  of  lth  panel  for  fuselage  (J=l)  or 

Jth  store  of  given  shape 

temporary  array;  contains  initial  velocity 
boundary  condition  which  is  destroyed  during 
the  solution 
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AN  { I , J )  temporary  array;  used  to  hold  the  aerodynamic 

influence  coefficients 

At  the  conclusion  of  the  calculation  of  the  fuselage  body 
source  strengths  identified  by  a  call  to  IOWRIT  in  Figure  A-2, 
arrays  XPT  through  GB  are  saved  on  TAPE7.  At  the  conclusion  of 
the  remaining  body  calculations  all  arrays  required  for  Program  II 
are  saved  on  TAPE10.  This  occurs  during  the  call  to  routine 
FRSTRT  identified  as  item  4  in  Figure  A-2.  For  the  fuselage, 
all  variables  and  arrays  in  common  blocks  BGEOM,  BOPTNS,  BINLET, 
BINSHK,  BSHOCK ,  DIMENS,  HEAD,  PARAM,  and  the  blank  common  arrays 
previously  saved  on  TAPE7  are  written  onto  TAPE10.  These  commons 
contain  all  the  information  necessary  to  restart  the  calculations 
in  Program  II. 

In  Program  I  when  both  a  noncircular  fuselage  and  an  elliptic 
store  are  to  be  analyzed,  the  results  of  the  fuselage  calculation 
are  saved  on  an  external  file  as  indicated  above.  The  data  for 
the  store  or  stores  overwrite  previous  values.  The  sequence 
employed  in  computing  elliptic  store  panel  strengths  is  performed 
in  a  slightly  different  manner  from  the  fuselage.  In  order  to 
take  maximum  advantage  of  similarities  in  store  shapes,  the 
geometric  properties  of  the  panels  and  aerodynamic  influence 
coefficients  of  a  given  shape  are  computed  only  once.  The  place¬ 
ment  of  a  store  in  a  different  position  or  orientation  causes  only 
a  variation  in  the  boundary  conditions.  Only  the  strengths  of  the 
source  panels  associated  with  the  additional  stores  of  identical 
shape  have  to  be  saved.  Provision  for  two  elliptic  store  shapes 
is  included  in  the  program. 

The  computation  sequence  for  the  stores  has  been  modified  to 
read  the  geometric  store  input  by  shape  first.  The  aerodynamic 
influence  calculations  are  then  made  for  each  individual  store. 

The  impact  of  this  sequence  of  computations  for  the  store  in 
blank  common  is  to  break  the  sequence  1  through  4  as  used  for 
the  fuselage  after  item  2  in  order  to  temporarily  save  all  arrays 
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on  TAPE11 .  This  break  in  the  sequence  is  indicated  by  item  5  in 
Figure  A-2.  At  this  point  arrays  XPT  through  ZC  and  common  blocks 
BGEOM,  BOPTNS,  BSHOCK,  DIMENS,  HEAD,  and  PARAM  are  saved  for  each 
of  the  elliptic  store  shapes  as  previously  described  for  the 
fuselage.  The  loop  on  the  number  of  shapes  is  then  repeated  in 
which  the  above  variables  are  read  and  the  source  strengths 
computed  for  each  of  the  stores  of  that  shape.  The  above  labeled 
common  blocks  and  arrays  XPT  through  GB  are  saved  on  TAPE10  by 
routine  FRSTRT.  Blank  common  is  then  used  to  copy  the  aerodynamic 
influence  coefficients  from  TAPE 9  and  TAPE8  to  TAPE10.  Blank 
common  is  equivalently  dimensioned  as 

COMMON  AN (KRAD, KRAD) , . . . , A (NADIM) 
for  the  TAPE 9  transfer  and  then 

COMMON  UB (3, KRAD, KRAD) , . . . , A (NADIM) 
for  the  TAPE 8  transfer  of  data. 

TAPE10  is  rewound  and  the  records  containing  the  fuselage 
arrays  are  read  for  use  by  DP HRS  at  item  6  in  Figure  A-2.  FRSTRT 
reinitializes  blank  common  and  the  various  labeled  commons 
identical  to  that  configuration  at  the  point  identified  by  item  3 
to  allow  use  of  the  noncircular  fuselage  solution  for  the  wing 
boundary  conditions. 

The  second  use  of  blank  common  in  Program  I  is  to  contain 
the  aerodynamic  influence  coefficients  associated  with  the  u-velocity 
panel  solution  for  the  wing-fuselage-pylon  of  the  parent  aircraft. 

At  item  7  in  Figure  A-2  DPCOEF  uses  blank  common  to  hold  the  coef¬ 
ficients  and  the  right  hand  side  through  the  panel  strength 
solution.  For  these  computations  blank  common  has  the  dimensions 
of 

COMMON  FVN (200, N+l ) 

where 
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FVN (I, J) 


u-velocity  panel  influence  coefficients 


FVN ( I , N+ 1 )  on  input  is  the  boundary  condition  of  the  Ith 

panel;  on  output  is  the  strength  of  the  Ith 
panel 

N  number  of  u-velocity  panels 

The  third  use  of  blank  common  in  Program  I  is  to  contain  the 
arrays  used  in  the  condensation  of  u-velocity  panel  strengths  to 
strengths  associated  with  single  corners.  Blank  common  is  used 
to  contain  the  arrays  of  corner  points,  net  strengths,  and 
orientation  information.  These  calculations  are  performed  in 
routine  NULYT  identified  as  item  8  in  Figure  A-2.  The  information 
is  retained  in  blank  common  only  until  it  can  be  written  onto 
TAPE12  by  WRFILE .  During  this  time  blank  common  has  the  dimensions 
of 


COMMON  XPT ( 500 ) ,YPT(100) ,ZPT(100) , SPHI (40) ,CPHI (40) , 

SWP (500) ,SNBP (40) ,CSBP (40) , THTBP (40) ,DPNET(500) , 
XPTS(1000) , YPTS (100) ,ZPTS(100) ,SPHS(40) ,CPHS(40) , 
SWPS (1000) , THTNET ( 1000 ) 

where 


XPT (I)  xw  coordinate  of  Ith  wing/body/pylon  constant 

u-velocity  corner  point 

YPT(J)  yw  coordinate  of  Jth  chordwise  row  associated 

with  wing/body/pylon  constant  u-velocity 
corner  points 

ZPT(J)  zw  coordinate  of  Jth  chordwise  row  associated 

with  wing/body/pylon  constant  u-velocity 
corner  points 

SPHI (J)  sine  of  dihedral  angle  associated  with  Jth  row 

of  wing  constant  u-velocity  corner  points 

CPHI(J)  cosine  of  dihedral  angle  associated  with  Jth 

row  of  wing  constant  u-velocity  corner  points 


SWP(J)  leading-edge  slope  of  semi-infinite  influencing 

triangle  associated  with  Jth  constant  u-velocity 
corner  point 
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SNBP (J) 

sine  of  orientation  angle  of  Jth  row  of  body 
interference  shell  u-velocity  corner  points 

CSBP (J) 

cosine  of  orientation  angle  of  Jth  row  of  body 
interference  shell  u-velocity  corner  points 

THTBP (J) 

polar  angle  in  cross-sectional  plane  defining 
the  Ith  row  of  fuselage  constant  u-velocity 
corner  points;  positive  in  counterclockwise 
rotation  from  positive  yB  axis 

DPNET ( I ) 

net  strength  of  Ith  constant  u-velocity  corner 
point 

XPTS  (I) 

xw  coordinate  of  Ith  wing/pylon  thickness 
source  panel  corner  point 

YPTS (J) 

yw  coordinate  of  Jth  chordwise  row  of  wing/pylon 
thickness  source  panel  corner  points 

ZPTS (J) 

zw  coordinate  of  Jth  chordwise  row  of  wing/pylon 
thickness  source  panel  corner  points 

SPHS (J) 

sine  of  dihedral  angle  associated  with  Jth  row 
of  wing  source  panel  corner  points 

CPHS(J) 

cosine  of  dihedral  angle  associated  with  Jth 
row  of  wing  source  panel  corner  points 

SWPS (I) 

leading-edge  slope  of  semi-infinite  influencing 
triangle  associated  with  Ith  source  panel 
corner  point 

THTNET (I) 

net  strength  of  Ith  thickness  source  panel 

corner  point 

Blank  common  is  last  used  in  Program  I  by  FRSTRT  as  previously 
mentioned  as  temporary  array  space  to  copy  the  source  panel  data 
for  the  fuselage  and  store  bodies  from  TAPElO  to  TAPE12.  See 
descriptions  of  FRSTRT  for  sequence  of  operations  involved  in 
file  transfer.  This  reference  to  blank  common  is  identified  as 
item  9  in  Figure  A- 2. 
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